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Definition 
drag coefficient 
lift coefficient 
moment coefficient 
tangential force coefficient 
moment coefficient due to normal force 

moment coefficient due to tangential force 
normal force coefficient 
pressure coefficient 

cavitation coefficient, base pressure coefficient 

constant part of the perturbation potential at the 
separation point 

series approximation of the Kernel function at the 
singular point 

nond imens ional normal component of the undisturbed 
velocity component (gn/U^) 
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nond imens ional strength of the source rings (q n /U oo ) 

velocity of the undisturbed flow field 
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angle between contour tangent and body axis 
nond imens ional potential function ( 0 / (U R) ; 
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body contourline from stagnation to separation point 
contour of the free streamline 
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radius of curvature 

density of the flowing medium 
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(n(s,a)/uj 

separation angle of the sphere measured from the 
stagnation point 
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SUMMARY 


Supercavitating flows about axisymmetr ic bodies at zero and small 
angles of attack are of wide practical importance. These flows are not 
so amenable to theoretical treatment as the plane flows, because the 
conformal mapping technique is not available in axisymmetric and three- 
dimensional flows. Approximate methods must therefore be adopted for 
the free streamline analysis. 

In this study an integral equation has been applied to calculate 
the free streamline flow behind axisymmetric bodies at zero and small 
angles of attack. The problem could be reduced to the solution of the 
mixed boundary value problem of potential theory with the additional 
condition that a specified velocity shall be attained along the free 
s treaml ine . 


I. INTRODUCTION 


The flow field about a body of revolution is usually determined by 
placing singularities along the axis of symmetry. Rankine [1] was thus 
able to calculate families of bodies by distributing point sources of 
variable strengths on the axis. This method was later improved by 
Taylor [2] and Fuhrmann [3]. Many years later, Weinstein [4] applied 
surface singularities to determine the flow fields about obstacles such 
as rings, disks, and cylinders. Further contributions to the indirect 
problem, in which the source strength is given and the corresponding 
body shape is sought, were made by Van Tuyl [5] and Sadowsky and 
Sternberg [6]. 

The direct problem, where the contour of the body is given and the 
corresponding strength of the sources is to be determined, was first 
solved by von Karman [7] with an axial distribution of singularities 
for the zero and finite angle-of-attack case. Lotz [8] probably 
published the first method using surface singularities; however, before 
that, Trefftz [9] had used annular sources to calculate the contraction 
coefficient of the vena contracta. Later, Riegels [10] extended the 
method of reference 8 and applied it to bodies of revolution that 


deviate slightly from rotational symmetry. All of these methods dealt 
with Dirichlet flows that do not permit the calculation of any drag 
force. 

The Helmholtz flow concept, since it enabled the aerodynamicis t to 
calculate a drag coefficient, meant an improvement. Over the years, a 
tremendous amount of literature on the two-dimensional Helmholtz- flow 
has accumulated. The first approximate solution of the axisymmetric 
problem was published by Bauer [11], who applied an axial distribution 
of sources to determine the drag of a sphere. However, this method is 
not well suited, since the flow field is everywhere continuous except 
on the a^cis . Therefore, the proper discontinuities cannot be duplicated 
at the separation point of the flow. 

Consequently, Armstrong and Dunham [12] applied a surface distribu- 
tion and devised an iterative scheme to determine the proper location of 
the free streamline. However, results calculated with this method were 
never published. The first numerical results for a disk placed normally 
to the stream were given by Garabedian [13]. 

In this report, the method of Riegels [10] is extended to the 
Helmhol tz-f low concept. The surface of the forebody and the free 
streamline are replaced by surface singularities. The location of the 
free streamline must be assumed for the first iteration. The final 
proper location must be obtained by an iteration that satisfies certain 
specified boundary conditions. The problem contains at least one param- 
eter for a forebody with fixed separation, probleme du sillage, the base 
pressure Cpg right behind the separation point, which is assumed to be 
constant along the remainder of the streamline. However, if the flow 
separates from a smooth body, for example, a sphere or ellipsoid, the 
problem will have one additional parameter, the separation point Sp, 
probleme de la proue. In cavitational flow, the base pressure is 
approximately equal to the vapor pressure of the liquid, and the pres- 
sure is constant throughout the cavity. For wake flow, the base pressure 
must be obtained from experimental data, and the assumption that the 
pressure is constant along the wake is questionable. Boundary layer 
theory can be applied to determine the separation point. 

A new linearized model is developed for the lifting body. The 
assumption of the Helmhol tz-f low is also valid for this model. The 
pressure within the separated cavity remains constant, and additional 
pressures due to the normal flow are equalized across the circumference 
of the cavity, thus causing the local normal force along the cavity to 
vanish. The cavity is therefore shifted to an asymmetrical location 
to satisfy the no-lift condition. 
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1 1. THE D IFFERENT CAV ITY OR WAKE MODELS 


Several free streamline models discussed in this chapter will be 
compared with the original Helmhol tz-Kir chhof f classical model. The four 
models are all characterized by the fact that the tangential velocity Vq> 
along the free streamline is greater than the approach velocity U^. 

For the Helmhol tz-Kirchhoff model, or infinite cavity model, the 
base pressure Cp-g is zero, and consequently the calculated drag is less 
than the drag observed in experiments. To correct this deficiency, 
higher velocities must be permitted on the free streamline. Then, how- 
ever, the streamlines curve back to the axis of symmetry, and certain 
cavity closure devices must be introduced. Thus, in model 2, the image 
or Riabouchinsky model, the free streamline reattaches to an artificial 
image of the forebody introduced at the end of the cavity. Only the 
drag of the first body is determined. The ultimate wake thickness is 
zero. In the third model, the reentrant jet model, the wake or cavity 
ends in a free stagnation point from which a reentrant jet projects 
forward toward the body base and vanishes there. This is an unreal 
feature, though it has some similarity to the often forward - thrown 
spray observed in cavities at low pressure. As in the Riabouchinsky 
model, there is a stagnation point at the end of the cavity, and the 
ultimate wake thickness is slightly negative on account of the fluid 
removed in the reentrant jet. One model, sometimes called the dissipa- 
tion model or parallel streamline model, probably describes the wake 
flow better than the cavity flow. Here, the downstream wake thickness 
is not zero. The pressure is initially constant along the streamline, 
springing from the disk edges until they reach their maximum wake diam- 
eter. From here on, the direction of the velocity vector remains con- 
stant; its magnitude, however, decreases until V^(s) is equal to U^, 
the velocity of the undisturbed flow. This model was used for most of 
the calculated zero-angle-of -attack cases of this report. However, it 
was used exclusively for the lifting cases. 

One can devise another model with exactly the wake width that is 
required to produce the drag force of the forebody. Here the far wake 
consists, again, of a parallel stream tube of radius 2. At some 

location downstream of the base of the body, the free streamlines inter- 
sect the stream tube forming thus a free stagnation point; therefore, 
the constant pressure condition cannot be satisfied over the rearward 
part of the free streamline, and other assumptions have to be made there. 

The literature abounds with other models which are more or less of 
practical usefulness. Originally all models were devised for the plane 
case; however, one can apply them to the axisymmetric case without 
difficulty. 
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III. DISCUSS ION OF THE THEORY 


The Potential and the Velocities of an Axisymmetric Body 

The perturbation potential function of an axis ymme trie body 
covered with surface singularities of strength q(s,w) = cos (noj)q n (s ) 
is given by the surface integral of the source strength times the 
reciprocal distance between the fixed (x,r,w) and the running point 
(£, P, w' ) : 


*E 2jt 


cp n (s,u),a) = 


JL_ 

4jt 



, v p cos (W ) dco' da 

q n O,a) K i L — 

[ (x- i) 2 +r 2 +p 2 -2rp cos (co-w' ) l 1 ' 2 


( 1 . 1 ) 


For the case n = 0, the body is covered with source rings of con- 
stant circumferential strength, representing the zero angle of attack 
case, and for n = 1, the source strength varies with the cosine of 
the meridian angle a). This case represents the body placed with its 
axis normal to the stream. 

In general, the strength of the source rings and the surface of 
the body is defined as 


q n (s,co,a) = cos (nco) q(s,a). 


To integrate the potential function around the body, we rearrange equa- 
tion (1.1) and write 


cp n (s, 0 J 5 Q!) = - ^ J q n (a, a) 


2jr 


cos (nco* ) 


dco 1 da, 


\ / (x-£) 2 +(rH-p) 2 o N/l-k 2 (l+cos (co f -w))/2 


where k 2 is the elliptic modulus 


k 2 = 


4rp 

(x-l ) 2 + (r+p ) 2 • 


( 1 . 2 ) 
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With the substitution 


to* - co = 2X; doj 1 = 2dX, 
we express the second integral in (1.2) as 


2jt 

\ ” f 7i 


cos (nco 1 ) dco 1 


o 


n/ l-k 2 (l+cos [co 1 - co] ) / 2 


= 2 f CQS [n(2X+o )}± dx< 
^ *j 1-k 2 cos 2 X 


(1.3) 


(1.4) 


To reduce the elliptic integral to Legendre's form, we use the substi- 
tution 


cos X = sin cp, X = arc cos (sin cp) 



(1.5) 


Reversing the limits -jt/2 and rt/2 yields finally, for (1.4), 

I„ - 2 {cos<™> 2 f l2n(,/2- cpU d(p 

n ^ 'il-k 2 sin 2 cp 


it/ 2 

s in (nw) J' 
-at/2 


sin f2n(rt/2-cp) ] 
1 -k 2 s in 2 cp 



With the relations 


cos [2n(^ - cp)] = (-l) n cos(2ncp) 

and 


( 1 . 6 ) 


sin [2n(| - cp)] = -(-1)° sin(2ncp), 
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one can define two integrals. The first one is 


it/ 2 

F (k 2 ) = (-l) n f cos < 2 ° jg dtp. 

^ \/l-k 2 sin 2 cp 


The second integral 


n/2 

d9 . o 

v 1 -k 2 sin 2 cp 


vanishes for all n. The perturbation potential can finally be written 
as 



<p„ <•. f 2 f V l!) ‘' 


^(x-|) 2 + (r+p) 2 


(1.7) 


Partial differentiation of the perturbation potential (1.1) with 
respect to x yields the perturbation velocity in the x-direction. 


^n 

dx 


s E 


P(x-E) 

[(x-|) 2 + (r+p) 2 ] 3 / 2 


2n 



o 


cos (nto 1 ) dco 1 

[l-h 2 (l+cos (oj'-co))/2] 3 / s 


( 1 . 8 ) 

With the aid of the substitutions (1.3) and (1.5), we can write for 
the second integral on the right-hand side of (1.8) 


2n 


%x 


cos (nco 1 ) dco 1 


it/ 2 


[l-h- 2 (l+cos (co' -to))/ 2] 3 / 2 


= 2 cos (nco) (-1) 


n 


SOS+2U& ^ 

[l-k 2 sin 2 cp] 3 / 2 


(equation continued on next page) 
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Til 2 

+ sin(nw) (-l) n 

-Till 



sin(2ncp) dcp 

[1-k 2 sin 2 cp] 3 / 2 


We define now the first integral on the right-hand side as 


V kE) 


x/2 



o 


cos(2ncp) dcp 
[1-k 2 sin 2 cp] 3 / 2 


The second integral on the right-hand side 


(- 1 ) 


Til 2 


-Till 


sin (2ncp) 
[1-k 2 sin 2 cp] 3 ^ 2 


dcp = 0 


(1.9) 


( 1 . 10 ) 


vanishes for all n. * 

The perturbation velocity component in the x-direction is now 


dx 



q n (s,a) 2p(x- g) G n (k 2 ) 
[(x-|) 2 + (r+p) 2 ] 3 / 2 


Partial differentiation of equation (1.1) with respect to the radius r 
yields the perturbation velocity component in this direction. 


dcp 

dr 



q n (a,a) • ■ - ■ — — 

[(x-|) 2 + (r+p) 2 ] 3 / 2 


2n 


o 


cos (nco 1 ) - p cos (uf -uS) cos (noo* ) 
[1 - k 2 (l + cos (w 1 -co))/2] 3 / 2 


dco' da. 
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The second integral on the right-hand side, together with substitution 
(1.3), yields 




cos (2nX + noo) 
[1-k 2 cos 2 X] 3 / 2 


dX 


rc 

2 P f co s . (2 jQ_ cos_ C 2nX _ +_n ( 4 dx> (1<n) 
c' [1 - k 1 2 cos 2 X] 3 ^ 2 


The numerator of the second integral can be written as 


cos (2X) cos (2nX + noo) = cos (noo) [cos {2(n-l)x} + cos {2(n + 1 ) X } ] 

- sin (noo) [sin {2(n + 1)X} + sin (2(n-l)X)]. 

( 1 . 12 ) 

With the definition (1.9) and the substitution (1.3), the first integral 
on the right-hand side of (1.11) becomes 


2 cos (noo) G n (k 2 ) 



cos {2nX + noo} 
[1-k 2 cos 2 X] 3 ^ 2 


dX, 


and the second integral is 


' cos (2nX + noo) 
(1-k 2 cos 2 cp) 3 / 2 


cos (2X) dX 


cos (noo) P cos [2(n-l)X]+cos [2(n+l)X] ^ 


(1 4- k 2 cos 2 X) 


3/2 


s in (noo) 
2 


1 sin [2 (n+l)X]+sin[2(n-l)X] ^ 

[ 14 k 2 cos 2 X] 2 


n/2 


__ 2 cos (no 0 ) 
2 


(-l) n 1 cos [ 2 (n-l)cp] + (-l) n+1 cos [ 2 (n-H)cp] 


[1-k 2 sin 2 ' 


: cp] 3 / ; 


dcp 


(continued on 
next page) 
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= 2 cos (nco) 


G (k 2 ) + G (k 2 ) 
n- n+x 


The integral with the factor sin (nco) vanishes for all n according to 
definition (1.10). The perturbation velocity component finally becomes 


^Pn = 1 

Sr 2jt 




= T- / q (a, a) 


2p cos (nco) 


[(x-l) 2 + (r+p) 2 ] 3 / 2 


f G n 


(k 2 ) 


_Q 

2 


[ G n . i ( k2 )4G ji+ i(k= 


)]ja< 


The perturbation velocity component in circumferential direction co is 
obtained by partial differentiation with respect to cj and multiplication 
with the factor 1/r: 


Sep 

Sw 


4jt 


5 2 tc 

d 2 P cos ( no)') sin(aj'-oj) 

q (a, a) ■ y / . " . 

11 [x-|) 2 +(r+p) 2 ] 3 / 2 ~ [l-k 2 (l+cos (co' -aj))/2] 3 / 


dco' da. 


(1.13) 


We apply (1.3) to the second integral in (1.13) and write 



2n 



cos (nw' ) sin(qj'-tjj) dco 1 
[l-k 2 (l+cos(co' -w))/2] 3 / 2 


1 cos (2atX+npj) sin 2X ^ 
[l-k 2 cos 2 X] 3 / 2 


(1.14) 


Using substitution (1.15) in equation (1.14) produces 

it/ 2 _ , i 

-2 sin(ncj) P (-1) cos (2 fn-llcp)- (-1) cos (2 (n+l)cp) ^ 
^ J (l-k 2 sin 2 cp) 3 / 2 

= sin (nco) [G n _ l (k 2 ) - G n+1 (k 2 )]. 
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The integral with the factor cos (nco) vanishes according to (1.10) and 
finally the w-per turbation velocity component becomes 


= J_ 

rSu) 2jt 



(a, a) 


o 


P 2 sin (nqj) 

[(x-D 2 + O+p ) 2 ] 3 / 2 



For zero angle of attack with n = 0, we obtain from definition (1.9) 
the relation G_ 1 (k 2 ) = G +1 (k 2 ). Therefore, dcp n /r^w equals zero for this 
case. The complete elliptic integrals F n (k 2 ) and G n (k 2 ) will be developed 
in appendix A. 

The derived perturbation velocities u, v, and w are the velocities 
for fixed points (x,r,oj) which do not coincide with the surface on which 
the sources are placed. The general solution of the problem, however, 
requires that we know the normal or tangential velocities on the bound- 
aries of the outer flow field. Part of the boundary is provided by the 
surface of the body. A closer investigation of the normal component of 
the velocity on the body surface reveals a discontinuity [ 141 . The 
jump in the velocity across the surface is exactly equal to the source 
strength, -hq(s,a)/2 to the outer flow field 5 and -q(s,a)/2 to the inner 
flow field. The tangential velocity component, on the other hand, is 
continuous. In this respect, we have to account for this discontinuity 
by adding the proper component to the various perturbation velocities 
on the body surface. 

We obtain for the axisymmetr ic flow, n = 0, 


q (s)r r (s) 1 p 2p(x-£) G (k 2 ) 

u (s) = - h / q (a) 

° 2 ° [(x-|) 2 + (r+p) 2 ] 3 / 2 


da, (1.15) 


and 


q„(s)x'(s) i P 2p[rG (k 2 ) - PG x (k 2 )] 

v (s) = / q (a) da. 

° 2 o [(x-i) 2 + (r+p) 2 ] 3 / 2 


(1.16) 
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The perturbation velocity components for the angle of attack case (n = 1) 
are 


u x (s,a) q-L^oOr' 1 p 2p(x-|) G x (k 2 ) 


OL cos go 


hf q ^ a - 


a) 


[(x-6) 2 + (r+p) 2 ] 3 / 2 


da, (1.17) 


v^s.a) q^s.aOx' 1 

n 


(X cos oj 




a) 


2p[rG 1 (k 2 ) - | (G 0 (k 2 ) + G 2 (k 2 ))] 
[(x-i) 2 + (r+p) 2 ] 3 / 2 


da 


(1.18) 


and 


w 1 (s,a) 

a(coS (jo) 


r, 


a) 


2p 2 [G Q (k :P ) - G 2 (k 2 )]/2 
[(x-6) 2 + (r+p) 2 ] 3 / 2 


(1.19) 


The notations r * (s) and x' (s) are differentiations with respect to the 
arc length s. For k~ = 1, the elliptic integrals have a pole with a 
logarithmic discontinuity. Therefore, the finite part of the integrals 
of equations (1.15) through (1.19) have to be taken. 

The total potential of the body is obtained by superimposing the 
perturbation potential on the potential of the undisturbed flow field. 


0 n ( s 5 w,a) = x(s) cos a + r(s) sin a cos w + cp^(s , go, a) . 

From the boundary condition 
30 (s,w,o:) 

3^ " V N = °> 


( 1 . 20 ) 


( 1 . 21 ) 
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we obtain the source strength on the surface of the body. We differentiate 
(1.20) with respect to the normal v 


d0 n (s,u),c>!) 

dv 


5x(s) 

5v 


cos a + 


sin a cos 


CO + 


Scp n (s,co,a) 


With the equations of appendix C, we can write now 


V N = 0 = -r f (s) cos a + x 1 (s) sin a cos go - u^(s r 1 (s) + v^(s, a) x 1 (s) 

( 1 . 22 ) 

where u n (s,a) and v n (s,a) are the perturbation velocities of equations 
(1.15) to (1.18). The first two terms on the right-hand side of the 
equation are the velocities of the undisturbed flow field. 

Let us consider first the axially symmetric flow field with n = 0, 
a = 0. The normal component of the parallel flow is independent of the 
meridian angle go: 


g 0 (s,a) 

The normal component sin 
gi(s,a,cd) = 


= r' (s) cos a ~ r' (s). (1.23) 

a of the parallel flow is proportional to cos go: 
-sin a cos go x 1 (s) « -x f (s) a cos go. (1.24) 


In the future, we will mostly use the linearized form of the boundary 
conditions. Only once in a while we will refer to the exact expressions. 

We obtain the integral equation for determining the source strength 
on the surface of the body by replacing in (1.22) the terms u n (s,QJ) and 
v n (s,c0 by their respective expressions (1.15) through (1.18): 


q n (s,a) = 2g n (s,0!) + -t J q n (a,a) K q (s , a) da, (1.25) 

o 
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where the Kernel function depends only on the geometry of the body 


K q (s, 0 ) = — — 

[ (x-|) 2 +(r+p) 2 ] 3 / 2 


[rx' -(x-g)r' ] G n ( k 2 ) 


' *£ < G n-!< k2 ) + W k *»} ' 


(1.26) 


The source strength q n ( s jQ0 is the only unknown of the integral equa- 
tion (1.25); all other functions are known. The elliptic integrals of 
(1.26) have a pole at k 2 = 1. The singularity, however, has logarithmic 
character, and the quadratic integrabil ity for solving (1.25) exists. 
Special procedures must be applied to solve the integral equation numer- 
ically. Equation (1.25), a linear integral equation of the second kind, 
shows the dependency of the source strength q on the angle of attack a 
and the shape of the contour of the body. If we set n = 0, we are con- 
sidering the axisymmetric flow field, and we replace q n (s,oi) by the 
linearized form of (1.23). The source strength is then not dependent 
on a. However, for the angle of attack case with n = 1, we use (1.24) 
and we notice that q 1 (s,a) =o;q 1 (s). 

Asymptotic Development of the Velocity Components 

The integrands of the integral expression for the velocity components 
are now developed into a power series of small e = a - s: 


p = r + r'e + r"e 2 /2 + ... 

§ = x + x'e + x"e 2 / 2 + ... , 


where primed values represent derivatives with respect to the arc length 
s. To solve the integrals numerically, the following procedure is used 
which shall be demonstrated as an example with the u-component. We 
repeat equation (1.15) for n = 0: 


u 

o 


(s) 


q Q (s)r f 


Tn f% 


(a) 


2p 

[(x-l) 2 + (r+p) 2 ] 3 / 2 


(x-|) G Q (k 2 ) da. 
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The integrand has a singularity at x = However, the value of the 

integral is finite if the principal value of the integral is taken. 

To avoid the singularity, the series approximation of the integrand 
containing the constant term free of e is subtracted from the integrand. 
Then the expression under the integral sign is zero at the singular 
point s = a, and the integral can be treated numerically. By designat- 
ing the approximating series as f(s,a) one obtains, for the u - component : 


u o (s) = 


qQ (s > r ' 1 _ 

2 + 2jt 


[q o (cr) K(s , a) 


f (s,a)] da 4- 


f(s,a) da. 


( 2 . 1 ) 

The first integral on the right-hand side is now free of any singularity, 
and can be solved numerically. The second integral contains the singu- 
larity and it can be treated analytically. The principal value of the 
integral must be used. This is the general procedure which is applied to 
all integrals as long as they are singular. 

Continuing now with the general development, one obtains, after some 
algebra, the series f(s,a) for the potential cp 


■f (s,a) = -q n (s,a) 


In 


8r 
I s-crl 


- 2n 


The integration of this function yields 
s 


E 


' Tk I £ J S >°) da = " 


q n ( s ,a) [- 


9 


2rt 


i 8r 

s_ In -] 1 

E s-a 


+ s In |1 - ~| + (l-2n)s 


+ 


where n = 0, 1, respectively. 

For the velocity components, the respective functions are 


f u (s,c) = -q^s.a) 


a - s 


x 1 r'+rx" 
2r 


■q ] [ 1 (s,Q:)x , + 
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r 


f v (s > a ) = -q n ( s ,a:) 


r' 1 , 8r 


o“ In 
a-s 2r 


,s-„ (1«°) ) + r ' 2 + rr " 




f (s,a) = — q (s,a;) 
w r n 


In 


8r 


s-a 


- 2J + ... , 


and integration again yields 


hf‘ u <s - 


g) da = - 


q n (s,oi) r 

2^r 


1 1 1 1 E 1 . x r ' + rx" 

x'ln 1 + - 

's' 2 r E 


q^C^aO 


2k X ' S E + 




a) da = - 


q (s,a:) 


n 


2k 


r 1 In 1 - — + 


S E, , r' 2 + rr" 


2r 


■ 77 - (s In -| — r + s In |l - — | 2n s F 

2r \ E s - s 1 s 1 * 


\ “ 

-| - 2n s E J 


- q' (s,a)r' s /2 k + 
n n 


h /v- 


a) da = n 


q n (s,a) 

2 nr 


8r 

s In -| 

E s - s T 


4- s In 1 - — 


+ 


The function f (s , a) and the respective integral for determining the 
strength of the source distribution (1.25) are 


f q (s,G) = q n (s,a) 


In 


8r 


s-a 


- (l+2n) 


r/c 


/2r + 
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Ei 

i/ 


f q (s,or) da 


q n (s,o:) 

2itr 


S E ln 


8r 


s-s T 


+ s In 1 - — - 2n s. 


- rK S E 


] 


+ • • . 


y 


where K is the curvature of the contour or the free streamline 


K - x 1 r ,f - r f x 11 . 


For acute cone apex angles, the above series expansion is not very well 
suited. For this case, a special expansion which is only applicable to 
small cone angles is given in appendix B. 

With the above series expansion, we can finally write the complete 
perturbation potential for arbitrary meridian angles go 


<P n (s,w,a) -I PC 2p F n (k2) / g r 

t — r — = - -r- / -s q (a. a) -- ■ - q (s,a) ( In - I I - 2 n })- da 

cos (n go) 2k J rn' r~ ~ . — , ttt n n v V I s-al 


x-|) 2 + (r+pp" u 


>} 


q n (s,a) 

2« 


n 8r 

ln -i 

E s-s T 


+ s ln 1 - — + (l-2n) s £ 


(2.2) 


The velocity components become 


u n (s,w,a) q n (s,a)r' 1 


cos (noo) 


+ 2 n J 1 q n (CT ’ a) 


2P(x-|) G n (k 2 ) 


[(x-£) 2 +(r+p) 2 ] 3 / 2 


+ q n (s,a) 


x' x 1 r 1 +rx n 

a-s 2r 


+ q 1 (s ,a)x f V da 


q n (s,a) 

~Tk 


x' In 1 - — + 

s 


Ei . x'r'+rx" 


2r E 


q^(s,a) 


2 it 


x V 


(2.3) 
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v n (s,w,a) q n (s,a) 

cos (nu) 2 


x 1 + 


u. 

A? / 


2p 


[(x-|) 2 +(r+p) 2 ]3/2 


( rG n 


(k 2 ) 


" 2 [G n-i (k2) + G n+i (k ^ ] ) + V 8 * 


a) 


- _^( ln 8r_ 
a-s 2r V a-s 


- (l+2n) 


) 


r'^+rr" ~| . t , . .1 9 °^ f . s t? . 

— J + q;(s,a)r'|da - ■ 2jt r’ In |l - -f\ + r 2 + rr 


-h(: 


S in 8r_ 
E s-s 


+ s In 1 - -S. 

1 s 


- 2n s 


.)] 


q' n (s,a) 

2 jc 


V 


(2.4) 


and 


w(s , U,(X) 
s in (nw) 


E 


n / v ( -> 8r 

- 7 q n (s,a) ( In 


CJ-S 


- 2 )} 


da + n 


q n (s,a) 

~ Sr 


8r 


S-r, In | , 

E s-a 


+ s In 1 


- 4] • 


(2.5) 


The integral equation (1.25) becomes now 


2g n (s,w,a) 


S n V B »w,ui; iff 

cos (nto) q n^ S,Q ^ + 7 J 


2p 


(a, a) 

o [(x-|) 2 +(r+p) 2 ] 3 / 2 L 


(rx'-r' (x-|))G n (k 2 ) 


(equation continued on next page) 
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px' 


_ (G (k 2 ) + G (k 2 )) 
2 n-i y n+i " 


q n ( s ,a) 

2r 


x' (In 


8r 


s-a 


- (l+2n) 


- r/c 


da + 


q n (s,a) 

2jtr 


g 

x* ( s In n — — — r + s In 1 1 - —I - 2n 

E s ~ s__ 1 s 1 E 


KX s T 


( 2 . 6 ) 


In this report, we very often refer to the normal and tangential veloc- 
ities of the body surface, which are quoted here for n = 0 and n = 1 
separately. According to the equations (C.l) and (C.2) of appendix C, 
we obtain for the normal velocity of the body in axial flow 


V s) ■ 


-r ' + 


q~( s ) 


2 - + 2 n J 1 2p V G > 


(rx'-r' (x-5))G o (k 2 )-x' pGiCk 2 ) 
[(x-|) 2 +(r+p) 2 ] 3 / 2 


V s > 

2r 


x 1 In 


8r 


cr-s 


x 1 - r/c ) y do 


+ 


%< s) 

2« 


fr v E ln 


8r 


s-s T 


+ s In 1 - — 


Ks , 


(2.7) 


The tangential velocity for this case is given as 
s 

p r (x' (x-|)+rr')G (k 2 )-pr'G 1 (k 2 ) 

V (s) = x' + — / i2pq (a) — — 

T 2rt J L ° [ (x- |) 2 +(r+p) 2 ] 3 / 2 


+ q o (s) 


+ 2 - ln , 

a-s 2r V |s-a 


q Q (s) 

- q ;<s)j-d 0 - 


pt) ln |l - ^| + ^ (l - ln T-5L 
2r J 1 s' 2r V s-s T 


q' (s)s„ 
H o v _ E 

2 it 


( 2 . 8 ) 
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The corresponding velocities for the body at an angle of attack CC 
are, according to (C.17) and C.18) for the meridian angle go = 0, 


lb 

v B1 (»,a) -a*' + 1_ J 

o 

^Zpq^a ,a) [ (rx 1 -r' (x-|))G 1 (k 2 ) - x' p(G 0 (k 2 )+G 2 (k 2 ) ) / 2 ] 
[(x-|) 2 + (r+p) 2 ] 3/2 


q i (s ’ a) f , n 8r , , u , 

' x In i r - 3x' -r/c ) r do 


2r 


+ 


q-L^a) 


2 It 


a-s 


~~ ( S r In 
2r V E 


r-^-r + s In |1 - -£| - 2s w 
s-s^ 1 s' E 


K S T 


(2.9) 


and 


1 

v T1 (s ,a) = ar' + 2 ^ 2pq 1 (a,a) -- 



(x' (x-|)+rr')G 1 (k 2 )- p r' (G 0 (k 2 )+G a (k 2 ))/2 


[(x-£) 2 + (r+|) 2 ] 3/2 


+ qi(s,a:) 




+ 


qi(s,a) 


} 


da 


q.,(s,aO r/ 

■srv 


r 1 s 
2r 


In 


.^1 

s 1 


2r 



q'^s.aOs^, 

2n 


( 2 . 10 ) 
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The Solution of the Integral Equation 


We replace the integral equation (2.6) by a system of N-linear 
equations and solve this system for the N-unknown q^. From now on, 
we reserve the index v for the fixed point s on the contour and the 
free streamline (which is here considered a fixed surface, too) and 
the index p. for the running point a. For the numerical integration 
scheme, we use the quadrature procedure of Gauss, where the value of 
the function at certain prescribed abscissas s^ is multiplied with a 
certain weight A . We call = Al^a^ the modified weight. It is 
obtained by multiplying the original weight a^, which is normalized to 
1, with the sectional length AI^ as shown in figure 1. The total 
length of the body is therefore given as 

N M P 

A = / AL 
H La 1 

p=l m=l [i= 1 

where the number of body sections is designated by M and P is the number 
of points in one segment. This number does not change. It is the same 
for all segments. The total number of points is therefore given by 
N = M • P. 




After replacing in equation (2.6) the integral sign with the sum- 
mation sign and collecting all terms multiplied with q v , we now obtain 


2 <*n>v 


cos (nw) 


1 + 


2r yrt r v 


8r % 

r + s,,ln |l - — I - 2ns 

IV S E I V S v' E 


N , 


“We 


x' In 


8r 


P- : 


=1 


s -s 
1 V |_l' 


(l+2n) - r K 

' v v 


i (r x'-r' (x -x ))G (k 2 )-r x' (G (k 2 )+G , (k 2 ))/2 

1 \ f •) 2r v v v v n n vt-i U v n- 1 vp. n+i vja A 

n ' (! |-l r„ .. \ 2 , /„ n2t3/2 




= 1 


[x -x ) 2 + (r +r ) 2 ] 3 / 2 

1 V v V |j/ J 


(3.1) 
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The primed summation sign means that the term jj. = v is omitted in the 
summation process, since f Vj ^was originally defined as the series 
expansion of the Kernel function K y ^ at the singular point. The 
inclusion of the constant term in tfris expansion insures that the two 
terms under the summation signs cancel each other at the point p, = v. 

In equation (3.1) the source strength q v is the unknown. The left- 
hand side of the equation is known and is substituted by either equa- 
tions (1.23) or (1.24) depending on n = 0 or n = 1. In matrix notation 
we write the N-linear equation as 


(3.2) 


The elements of the main-diagonal are given by the expression 


a ll a 12 •• a lN 


q ni 


2 ( g n^ l 

a 21 a 22 • • a £|^ 


q n2 


2 (g n ) 2 

a, a >T . . 

Ni N2 NN 


q nN 


2 ^n>N 

_ 


_ 




a W 1 + 2jtr 


V *- 


8r v S E 

s -sj + s v ln I 1 - —I - 2ns ] 

V t 


Vv S E 


N . 


2 nr 


8r , 


x' ln 

v V s -s 


- (l+2n) - kv 


V V 


p.=l 


and the elements of the neighboring diagonals are 


a = 2r - — 

vn n 


( r v x v" r v ( G n-i (k v 2 n )+G n +1 ( k vu ))/2 A 


n[(x -x ) 2 + (r +r ) 2 ] 3 / 2 
v \1 V (J. 


U* 
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The matrix is usually conditioned satisfactorily [ 15 ]. The elements 
of the principal diagonal a vv are larger by one magnitude than the 
elements of the neighboring diagonals. The solution function q v 

is usually very smooth. For one particular example of calculation, 
the value of the determinant was det A = 53.4. The conditioning 
number of the system had a value of = 2.7 x 10~ 14 . The number of 
points used to calculate the flow field was N = 64. The determinant, 
as well as the conditioning number, depends on N. With decreasing N, 
the conditioning of the system usually improves, provided the deter- 
minant is different from zero. 

After having solved the system of equations for the source strength 
q v , one obtains the velocity components by a simple quadrature: 


(u n>v 
cos (nco) 


(qj, 


n , 


r' + 


2 v 2 tc 


i=1 




2r v (x -x )G (k 2 ) 

v v u n v y u / 

[(x -x ) 2 +(r +r ) 2 ] 3 / 2 
v p v p 


x 1 x 1 r ! +r x M ; 

4- (q ) V — + V V _ V V V+ (q ! ) x 1 

V4 n 7 v 1 s -s 2r J VM n v v 

p v v 




< q n>v 

2n 


s_ x r r'-fr x M 

X* In |1 - -^| + v v 9 v v s p 
v 1 s 1 2r E 

v v 


«i>v 


2« v E 


x'a , (3.3) 


(v ) 
v n v 

cos (nco) 


N i 


(q n ) v , . 1 \ I t v n 

2 X v + Tjt / , 1 (q nV 


i=l 


2r [r G (k 2 )-r (G (k 2 )+G v (k 2 ,))/2] 
p L v n v vp 7 p v n-i v vp 7 n+i v V P' 7 


[ (x -x ) 2 + (r +r ) 2 ] 3 / 2 
v p v p 


+ (q ) 
vn n y v 


s -s 2r 

L- P V V 


If V 

' In tt — r - (l+2n) + 


s -s 
P v 1 


r ! +r r” “• 
v v v 

2r 


(qj 


s r 1 +r r * * 

W f [ A | , , ,, E j v v v 

+ r T > A : ir In 1 1 - — | + 


2jt *■ vj "|J. 2it v 
8r, 


2r v E 


, s In -r \ + s In 1 1 -| - 2ns - r' s 

2r v V E IVeT v 1 s v « VJ 2 * vE 


(3.4) 
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and 


<Vv 

s in (nco) 


, r r 2 (G (kf )-G , (k 2 )) 

±- \ J(q ) n ~ 1 V(j, _ 

2n Aj l n V- [(x -x ) 2 +(r +r ) 2 ] 3 / 2 
u=l LV v w ' v V (/ 


^ <«»>v 


8r v Tl 

In r -2 /-A 

S -S 1(0. 

1 JJL V ^ 


(3.5) 


+ n 


2jtr,, 


8r v , S E, 

s In r + s In 1 — - s 

- E VV V E - 


This is the usual procedure to solve the flow field about a body of 
revolution with a fixed contour. So far, the wake or cavitation bubble 
has not been mentioned specifically. However, for the solution of the 
wake problem, which will be discussed in the following chapter, a solu- 
tion q v of (3.2) is necessary and will be used to calculate the perturba- 
tion potential and the tangential velocity. 

The Axisymmetric Cavity 

Following Helmholtz’s suggestion, the free streamline discontinuity 
surface is taken as an idealization of the separation surface which 
divides the main flow from the wake or cavity which follows separation. 
The thickness of the sheet is taken vanishingly small, and the flow on 
one side is assumed to have no effect on the other except through the 
pressure. For cavitating flows with ratios of cavity density vanish- 
ingly small, the flow in the cavity is assumed as quiescent, and it is 
therefore a constant pressure region. 

The problem to find solutions to the discontinuous flow about 
bodies of revolution with a smooth contour has two parameters. These 
are the separation point sp and the cavity pressure coefficient Cpg. In 
order to explain the general method of solution in principle, we cal- 
culate the flow about a body where the flow separates abruptly at the 
corner of the base. The problem then has only one parameter which is the 
cavity pressure coefficient Cpg or the equivalent tangential velocity V.p^ 
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These prescribed parameters are constant along the free streamline 
in accordance with the assumption that the flow in the cavity is 
quiescent. For small velocities in the cavity, Cpg and Vpg are func- 
tions of the arc length s. For the discussion of the general theory, 
we assume a constant Vq*g and we adopt the infinite cavity concept, 
which is the Helmholtz model (Vmg = 1), Special changes of the general 
method of solution if other models are applied will be described at the 
end of this chapter. 

In particular, we assume the flow to be (1) inviscid, (2) gravity- 
less, and (3) steady. Figure la represents the upper half-plane of 
a body of revolution and the corresponding cavity. According to our 
assumptions, the body possesses a sharp corner at sg where the flow 
separates abruptly. Otherwise, the contour r 2 °f the body is smooth 
except at the stagnation point. The cavity is separated from the 
outer flow field by the free streamline r^. On^the given contour 
Ti, the absolute value of the velocity vector |Vm(s)| is unknown; its 
direction, however, is known with condition (1.21). Downstream of Sg 
the absolute value of the tangential velocity vector | VxB I specified 
on r 2 ; its direction, however, is unknown. With the specified velocity 
Vtb on the velocity potential is also known on r 2 except for a 

constant. The position of r 2 is therefore not known and we have to 
choose it arbitrarily on some reasonable basis (see Fig. lb). We 
designate this line r*. 

It is our goal to determine the position of r 2 . For the solution 
of the problem, we formulate that: 

(a) The boundary condition Vn(s) = 0 has to be satisfied on 
Fi- 

(b) On the arbitrary line r*, the velocity and therefore 
also the potential except for a constant is known. 

(c) If r§ is not identical with the free streamline r 2 , 
the normal velocity (V^(s) £ 0) does not vanish on r 2 
and therefore is not a streamline. 

(d) From the magnitude of Vjq(s) on r% 9 we can deduce the 

shift of the line The solution of the problem is 

therefore obtained by iteration. 

(e) Since the position of the free streamline is not known, 
we have to conduct all numerical calculations on the 
arbitrary line 
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II 


Besides the velocity Vjg, we mentioned also the potential 0 ob( s ) 
in (b) . Later on, we will employ this potential function rather than 
the velocity V-pg. The reason for this choice will be given later in 
this report. 

(1) The Solution of the Problem by Specifying the 
Potential Function 

As a next step, we establish the set of integral equations with 
which we determine the source strength q Q (s). On the contour Fi, we 
satisfy' V^(s) = 0 and obtain therefore 

S E 

q (s) = 2r' (s) + \ f q (a) K (s,a) da, (4.1) 

o rt J o q 

o 


where the Kernel function Kq(s,a) is given for n = 0 by equation (1.26). 
For points within the region we obtain the integral equation 


0 aB (s ) - x'(s) = 


hi 


q (a) K (s,a) dcr 
o cp 


(4.2) 


with the Kernel function 


K (s,a) 
cp 


2pF o (k 2 ) 

I (x*- |) 2 +(r*+p) 2 


With known potential 0 oB (s), equation (4.2) is an integral equation of 
the first kind. In both equations (4.1) and (4.2), the source strength 
q (s) is the only unknown function. The Kernel functions Kq(s,a) and 
K^(s,a) depend only on the geometry of the body contour Ti and the 
arbitrary line r^. The integration limits are the stagnation point 
(s = 0) and the point where the wake is truncated. About the choice of 
the point s E and the influence of the truncated part of the "infinite" 
cavity or wake, some explanation will be given at a later time. 
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On the left-hand side of equation (4.2), we recognize the 
total potential of the free streamline p 2 which is defined as 


>V S > 


<*oB< S B> + f 


V™ dS > 

TB ’ 


B 


(4.3) 


where only the value of the integral is known. The value of the poten- 
tial at fehe separation point Sg is a constant and unknown. 


^oB^ S B^ 


C. 


The two terms on the left side of equation (4.2) would represent the 
perturbation potential of the body with the exact free streamline r 2 
if x*(s) would have been the abscissa of r 2 . Since, however, x*(s) 
belongs to the arbitrary line r %, the term 02f o ;g(s) - x* (s ) ) is only 
approximately cp 0 g(s), except at the separation point where x' v (s) = x(s) . 
For a convergent procedure, however, where r 2 = r 25 the left-hand side 
becomes, after a sufficient number of iterations, identical with cp oB (s). 
We define there 


cp oB ( s ) = i^oB ^ ' x ( s ) = <P 0 ( S ) + c - 


(4.4) 


The barred expression 


cp o ( s ) = J V TB ds " x “( s ) 

S B 


(4.5) 


represents all known functions. The subscript B refers here to quantities 
compatible with the specified pressure coefficient Cpg. 

The perturbation potential cp oB (s) can also be represented by 
the function cp*(s) plus some corrective term which will approach 0 as 
p* p 2 . We obtain cp*(s) by solving the system (3.2) for q*(s) and 
insertion of this function into equation (2.2), where we set n = 0, for 
the axisymmetric flow case. Since the system (3.2) was obtained with the 
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boundary condition V^j(s) = 0 along r x and V% > the line is temporarily 
a streamline. All functions obtained with q*(s) of the system (3.2) 
will in the future be marked by asterisks, in order to distinguish them 
from those obtained with the solution q Q (s) of the system of integral 
equations (4.1) and (4.2). We add and subtract in equation (4.3) the 
terms 


cs 

J2£(b b ) and J v*(s) ds 


and obtain 






0O <V + ' f V T (s) ds + J 


+ I ( V tb- v t ( s )) ds 


Due to the definition (1.20) of 0*(s), the terms 




Vj(s) 


ds = 


0>) 


= X x (s) + cpo(s) 


B 


can be collected, and we obtain instead of equation (4.4) now 


W* 0 = ^ (S) ‘ 


x*(s) 


/ 


(v tb" v t ( s)) ds + C = 9 0 < s > + c. 


(4.6) 


27 



where the constant is given as 


C ^oB^ S B^ ” ^oB^ S B^* (^- 7 ) 


The barred function in equation (4.4) can be expressed with cp' c (s) = 
0o( s ) ~ x*(s) as ° 


cp o ( s ) 




V^,(s)) ds. 


(4.8) 


All terms on the right-hand side are known quantities. We insert equa- 
tion (4.6) into the integral equation (4.2). Since cp Q g(s) is not neces- 
sarily identical with cp*(s), the normal velocities on are not zero 
(Vn(s) £ 0) , and r* is not a streamline any more. Our goal is to shift 
1% in such a fashion that (a) the condition of zero normal velocity along 
F* is satisfied (r% -> 1^) and (h) the specified perturbation potential 
cp, 0 B( s ) is attained along r*. 

The representation of the exact perturbation potential as a 
correction of the approximate cpo( s ) has certain advantages for the 
numerical calculation. The advantages will be discussed at the end 
of this chaper in more detail. 

The problem still has one unknown quantity: the constant C, 

as given by equation (4.7). We must therefore later establish an equa- 
tion for calculating C; only then the problem is completely determined. 
The constant C is for the time being arbitrarily assumed to be unity. 
However, since the geometric functions of are introduced into the 
Kernel functions of the integral equations (4.1) and (4.2) the exact 
free streamline can only be obtained by an iteration procedure. 

The system of integral equations (4.1) and (4.2) will now be 
satisfied in N points along the contour and the line r*. For this 
purpose, we replace the integral sign by a summation sign and apply the 
Gaussian quadrature. For simplicity, we drop the asterisks for the 
geometric functions of the line r**, but we keep in mind that x(s), r(s), 
and their higher derivatives belong to the assumed line JT*. Occasionally, 
we will return to the use of asterisks if it seems necessary to the 
understanding of the problem. 
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For the points 1 I v I i on the contour r 2 > we obtain the equation 


, __ v r 


r = 
v 


2tct 


N . 


8r 


i=l 


In -r — 

v V s -s 
1 n v 


r k >■ A 
v vj n 


*=: { x v (■ 


8r 


s t, 1° 1 

E s -s„ 
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+ s v In ll —I -tks 

v 1 s' v v E 


, ' 2 (q ) r 

. _L \ 

+ 2n L 


i=l 


2 (q ) r [(r x'-r'(x -x ))G 0 (k 2 )-r x'G x (k 2 )] 
n " " y v v v u ° vu u v x VLL 

[(x -X,,) 2 + (r +r ) 2 ] 3 / 2 


A . (4.9) 

M- 


v )i- 


V |1 


The point i is the last point on the body contour immediately upstream 
of the separation point sg. The numbering system of the points v is 
explained in Figure 1. 

For the points j ^ y < N on the line r'g, we obtain the equation: 




+ 1 = 


, \ N 

tao> v r \ 


2n 


8r 


In 


i=l 


s -s 

1 |1 V 1 


A - 


8r 


s„ In H 

E s -s 

v E 


+ s In 
v 



, * 2r F (k )A 

-h) M w H r • 

2rt ‘-i ° 11 [(x -x ) 2 + (r +r J 2 ] 1 / 2 

\ 1 = 1 v p/ v V p/ 


(4.10) 

where j is the point on r* immediately downstream of the separation 
point Sg. 


Equations (4.9) and (4.10) replace now equations (4.1) and 
(4.2), and form a system of N linear algebraic equations whose solu- 
tion is the source strength q Q (s) along and The equations are 

in matrix form: 
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a V, 
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ii 
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12 
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q - +Aq • 
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q - + Aq • 

OJ ^OJ 


2r 1 + 0 


2r. + 0 
1 


cp . + 1 

^°3 


( 4 . 11 ) 


q oN + A< *oN| J^oN + 1 ' 

The coefficients of the matrix are calculated by the expres- 


s ions ; 
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for the region j g v S N and 
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A r F (k 2 ) 

u u o vu. 

* [(x -X ) 2 + (r +r ) 2 ]i/2 

V [i V [I 

for the region j ^ v ^ N and v \i. 

An estimation of the magnitude of the elements of the matrix 
shows that the elements increase in size toward the main diagonal. 

Also, the sign of the elements is the same on either side of the main 
diagonal. The elements a v ^ on and T% are essentially represented 
by ln(8r v | Sy-SjjJ ) . We remember that the integrals possess a logarithmic 
singularity. The elements of the main diagonal are determined by the 
finite part of the singular integral. Ifi calculated examples, the main 
diagonal elements were greater than the neighboring a v ^ by at least a 
magnitude . 


a 

VU 


The right sides of the system of equations are formed now. 
According to equations (4.1), (4.2) and (4.4), we obtain 


and 


b = 2r* for the region 1 ^ v ^ i 
v v 


h v = 9, 


ov 


for the region j ^ y ^ N. 


The solution of the system with these right sides yields q QV . 
The solution q Q (s) usually jumps at the separation point. However, if 
q Q (s) is also the solution of the system (3.2), then the source strength 
q Q (s) behaves like r* (s) at this point; r* (s) has in s-g a corner for the 
abrupt and a contact point of first order for the smooth separation. 

If the source strength q 0 ~ is inserted into the equation for 
the tangential velocity, Vt v on is not necessarily equal to V'pg. We 
have to add an additional source strength Aq ov , which is multiplied with 
a constant C to obtain V^b on r^. We do not know yet the constant C of 
the perturbation potential; therefore, we set in (4.11) C = 1 along the 
line and on the body contour Tj, the constant is C = 0, since the 

normal velocities shall still vanish there. 

The right sides of the system (4.11) are now, with the same 
elements a y given as 


b^ = 0 for the region 1 ^ y ^ i 
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and 


b =* 1 for the region j ^ v ^ N. 

We solve this system for the additional source strength Z;q ov . This 
source strength refers to the constant C = 1 and always possesses a 
jump at s s , since the right sides have a jump there. 

For an arbitrary constant C, we obtain the complete solution 
of the system of equations (4.11), therefore, as 


o v 


- q + CZ\q . 
o v O v 


(4.12) 


Before we discuss the determination of the constant C, we want to list 
the quadratures for the tangential and normal velocities. Gauss quad- 
rature procedure is applied to equations (2.7) and (2.8), and we obtain 
two expressions which contain the derivative of the source strength 
^ov = (dq o OO/ds) an< 3 which can be combined: 
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The tangential velocity is thus given as 
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(4.14) 
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The normal velocity is given with equation (2.7) as 
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(4.15) 


We now turn to the determination of the constant C of equations (4.4) 
or (4.12). 


We prescribe on r| in v = j the tangential velocity to be 
equal to the given Vtb one. Since we have already corrected the 
perturbation potential of r* to that of the exact free streamline and 
obtained the corresponding q Q (s), we expect the tangential velocity 
along the entire length r ^ to be equal to V-pg. 

The complete solution (4.12) is inserted into equation (4.14), 
and we solve for the constant C. 
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where the function F x is given as 
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and the Kernel K VT (fj,j) is 




2r {(x*. (x.-x )+r.r')G (k 2 .).-r r'.G-^k 2 .)} 
?j) = _ J v 1 iL J-JT Q \kll2 H j lV 


[(x -x.) 2 + (r +r.) 2 ] 3 / 2 
V M. J [i J 


The differentiation of the source strength will be discussed later on 
in greater detail. The function £q 0 changes rapidly on r* in the 
neighborhood of Sg* It is therefore difficult to obtain the derivative 
(Z^q^)j sufficiently accurate. We encounter these difficulties always, 
during the first integration steps where C is relatively large. In 
this case, we use C as obtained by (4.16) as the zeroth approximation 
for determining C by the regula falsi with equation (4.14). 

After the constant C has been determined, we insert (4.12) 
into equation (4.15) and calculate the normal velocity component V^(s) 
on and improve the coordinates of the new line 

From the velocity diagram of figure lb, we obtain for small 
angular differences, the radial component of Vrp(s) as: 


r'**(s)V T (s) = r'*(s)V T (s) + x'*(s)V N (s) - r'*(s)v2(s)/ (2V T (s)) 

(4.17) 


and the axial component 


x'**(s)V T (s) = x'*(s)V T (s) + r'*(s)V N (s) 


x'*(s)V|(s)/(2V t (s)). 

(4.18) 


34 



Multiplication of equation (4.17) with x f *(s) and (4.18) with r ,#c (s) 
and subtracting the latter from the first yields 


V s > 

V T (s) 




(s) r 


'**(-> 


I * I 

r 1 x 




( 8 ). 


We solve the equations for x 1 ** and r 1 " and obtain the derivatives of 
the new line which will be much closer to the exact streamline f 2 
if the procedure converges. 


X'**(s) = x'*(s) 41 - (V N (s)/V T (s)) 2 - r'*(s) V N (s)/V T <s) 

(4.19) 

and 

r ' (s ) = r'*(s) 4l - (V N (s)/V T (s)) 2 + x' * (s) V N (s)/V T (s). 

(4.20) 

The second derivatives x r, ' f ' v (s) and r"'*(s) are obtained by a numerical 
differentiation of the corresponding slopes^ In the vicinity of the 
separation point s^, the slopes of r 2 ar *d r 2 " , respectively, are to be 
considered separately in order to find the proper shape there. Inte- 
grating the derivatives r f **(s) and x'*' v (s) yields the new coordinates 
of the improved line r 2 ' v 


r**(s) = r(s +e) 
B 



JI.JU 

' wf (s) ds 


and 


x**(s) - x(s^+e) + J' x ? *^(s) ds. 

S b +€ 

The determination of the coordinates in the region s^ ^ s ^ s^ + e will 
be discussed later. 
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Two possibilities of representing the perturbation potential 
ep 0 g( s ) were given at the beginning of this chapter. The reason why the 
perturbation potential (4.6) was chosen shall now be given: 

Ihe matrix of system (4.11) is in general not very well con- 
ditioned [l3]. For one specific case, which was the flow around the 
disk with Cpg = 0 and N = 64 points on Fi and r*, the determinant had 
the value of det A = 5.4 x 10“ 33 , and the conditioning number was 
Kg = 2.8 x 10~ 52 . If the right sides of the system were determined by 
equation (4.4) and r* in the region j ^ v ^ N (case A), then the func- 
tion (q Q ) v , as well as (A} 0 ) v , had a discontinuity of approximately 
equal magnitude at the separation points sg. The constant C was approxi- 
mately 1. Both functions had small oscillations in the vicinity of Sg 
on F% which did not disappear after the combination of both functions, 
according to (4.12). The amplitude of the oscillations was relatively 
small; therefore, the procedure still converged. However, the dif- 
ference in the tangential velocities (Vtb “ Vx(s)) did not fall below 
a tolerance of .01. These oscillations did not appreciably disappear 
when a larger number of significant digits were employed. 

However, if equation (4.6) determines the right side of the 
system (4.11) (case B) , the jump in the function (q Q ) v at Sg disappeared 
with equation (4.8). We remember that originally the perturbation func- 
tion (cpo) v was obtained with a continuous source strength (q*) v of the 
system (3.2); and (q*) v is continuous if r 1 (s) is continuous on and 
r%. In the second solution (Aq 0 ) v of the system, the oscillations 
remain, since the jump in the right side remains the same. However, 
the oscillations disappear now with the constant C. The tangential 
velocity difference (Vtb “ Vt(s)) could be made smaller than .001 for 
this case. 


Where accuracy is not the leading factor, approach (A) is 
the shortest to obtain a solution. Here, the additional solution of the 
system (3.2) is not called for, since cp*(s) an ^ V<x(s) do not appear in 
the evaluation of the perturbation potential. 


(2) The Direct Solution of the Problem by Specifying 

Vtb on 

Instead of using the specified potential along f'z, one under- 
standably would prefer to use the tangential velocity Vtb directly.. 

With this approach, one would eliminate the determination of the con- 
stant C, and the procedure would be shortened considerably. However, 
one has to overcome other obstacles. 
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In lieu of the 
now equation (4.14) with 


integral equation (4.2) or (4.10), 

(Vt)v = ( v TB)v for the region r*: 
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This expression represents formally an integro-dif fer ential equation for 
the unknown quantities (q Q ) and (q' Q ) v . The geometric functions x y> r yJ 
and their derivatives are temporarily provided by the line F*. We 
express the derivative of the source strength (q^) v by the differentiated 
Lagrange interpolation formula. 


(q*) = L* (q ) + L* (q ) + L f (q ) . . 

VH 0 V V-l CTy -1 V 0 / v V+l O V+l 


The coefficients of the interpolation formula L i (s) are polynomials of 
degree n, which contain only s^. 

The elements of the matrix (4.11) for points of the region 
(j ^ v ^ N) on r* are given now by 
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For the points v £ jj. and \± v + 1, the elements are given by 


a 

vja 


A r 


[x f (x -x )+r r f ]G (k 2 )-r r'G-tfk ) 
v v la v v o Vja p v v[x' 

[ (x -x ) 2 + (r +r ) 2 ] 3 / 2 
v |a y v \i 


For the diagonals immediately neighboring the principal 'diagonal , we 
write, with the above expression. 


+ 

a v(v+i) 


1j* - + a - . . 

v+i v ( v+i) 


If the points v = j and y = N are approached, the abscissas s v of 
Lagrange's interpolation formula have to be shifted downstream, or 
upstream by a point, respectively. In the neighborhood of the separa- 
tion point s-g, the differentiated interpolation polynomial represents 
only a rough approximation to (q q ) v , since (q^) v is singular there for 
the abrupt separation. 

The integrand of the integral expression for the tangential 
velocity component was developed into a series with small e = (s - Sg) . 
The dominating term was 1/(sjj-s v ). The elements a v ^ therefore change 
sign when passing across the main diagonal. The absolute value of the 
elements on both sides of the principal diagonal is large. In all cal- 
culated cases, the value of the elements |a v ( v -f :L )| was larger than the 
value of the element of the main diagonal, which was given by the finite 
part and the central differentiated coefficient of Lagrange's interpola- 
tion polynomial. A matrix of this kind is usually not very well condi- 
tioned. The resulting solution q Q (s) was oscillating severely; therefore 
it was not fit to be used for determining the normal velocity V^(s) with 
equation (4.15). The procedure diverged. The amplitude of the oscil- 
lation varied with the number of points along F*. Other methods to 
solve the system (4.11) were not employed. This method of solution was 



(3) The Free Streamline in the Neighborhood of the 
Separation Point 

The special behavior of the free streamline F 2 at the separa- 
tion point is now investigated in more detail. In analogy to the plane 
cavity flow, one can calculate families of convex and concave free 
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streamline flows, which have the separation point sg and the cavity 
pressure coefficient Cpg as parameters. One can also apply the 
behavior of the plane free streamline at the separation point to the 
rotational symmetric one, especially since the two first terms of the 
Laplace equation 0^ + 0 £r + 0^/v = 0 at and in the vicinity of a 
singularity are overwhelmingly large compared to the third term of 
the left side. In this respect the solutions of the Laplace equations 
for both flows are similar. This argument is definitely true for the 
abrupt separation. For this case the free streamline r 2 the infinite 
curvature at Sg and the tangential velocity has an infinite gradient 
on Tx at S B • 


Smooth bodies with continuous curvature, furthermore, contain 
points on their surface in which smooth separation occurs. For these 
points Sg the curvature k of the free streamline r 2 i s finite and equal 
to the curvature of the body at sg. The tangential velocity is smooth 
and continuous at Sg when passing from Fj to r 2 * The discussion of 
these flows will be taken up again at a later time when special cal- 
culated cases are presented. 

The free streamline is approximated by a suitable polynomial 
at the separation point sg. We assume, in accordance with the above 
argument, that the two-dimensional free streamline has the same 
behavior as the rotational one at Sg. Let us, therefore, for illus- 
tration, consider and discuss the Helmholtz flow around a plate of 
width b = 2. From reference 16, we obtain the coordinates of the 
free streamline r 2 as a function of the parameter t. 

x=^[t^T - In 1 1 + \/ t 2 -l | ] ; r = f^(t-l) + | . 

(4.21) 

The parameter t attains values between t = 1 and t = oo. For t = 1, we 
obtain the separation point s B of the plate. To change the equations 
(4.21) to the arc length s, we differentiate x and r with respect to 
t, form the expression ds/dt and integrate over t. The arc length is, 
consequently. 


s 



(t 2 -l). 


The coordinates finally become, with the abbreviation e = s-s B , 


r(s) 


1 
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and 


x(s) 


b 

4-fxt 





The first derivatives of the coordinates are 


x' (s) 




(4.22) 


and finally the second ones 


x" (s) 



r"(s) 




(4.23) 


These combine to yield the curvature 


/C(s) 


x" (s ) 
x'(s) 



Looking at the second derivative, we observe that the singularity is 
caused by x M (s). The expression r"(s), however, approaches with € -> 0 
a finite value. 


lim r" (s) = - “ = -1.7854. 

e-> 0 


For the rotational case, we choose a similar polynomial which, 
however, includes a wider variety of cases, the abrupt as well as the 
smooth separation. After having inspected a large number of cases, we 
decided that the expression 


r 1 + x* \/ £Sl x 

r 1 (s) = — — h r^e + b ± e 3 ^ 

n/ 1 + ea-j_ 


(4.24) 
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was the best fit to the free streamline r 2 in the vicinity of the separa- 
tion point S 33 . The terms r^, x^ and their higher derivatives are the 
slopes, etc., of the body contour at Sg, = x f (s-g) . 

According to the relation x f (s ) 2 = 1 - r f (s) 2 , we obtain, 
after some algebra. 
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x* ±r * Jea, 

o o 
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r 2 (r"e + b^ 3 ^ 2 ) (r *+x f \Tea^) 4 1-fea x 
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(4.25) 


The second derivatives are 
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+ x * \l a 7 a 3 (r f + x* sfeal) 
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(4.26) 


and, with r f r" -f x'x" = 0 , we obtain 


x M (s) - 


r 1 (s)r fl (s) 
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The second derivatives show again which of the expressions contribute to 
the singularity. In the case of a flow separating from the disk, with 
Xq = 0 and r^ = 1 , x ,! (s) causes the singularity, and r n (s) approaches a 
finite value (a-jr^/ 2 + r”) . If we set a x = (4+jt)/b, we obtain the exact 
expressions for the plane case (equation (4.23)). On the other hand, 
if x' Q = 1 and r q = 0, r"(s) will cause the singularity. This case, for 
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instance, represents the flow about a sphere with the point of separation 
at s B = it/2. 

We choose the sign according to 


r' (s) - (V + r'^e + b^e 3 / 2 ) J 0. 


We obtain this relation from equation (4.24) with the denominator of the 
first term neglected. At the point of smooth separation, the square 
root \f~e aj changes the sign and r' (s) has for this case the series 
development; 


r' + r"e + b-,e 3/2 + 
o o ± 


Integration of equation (4.24) yields the radius of the free 
streamline We integrate between the limits s-g and s, or e = 0 and 
€, and obtain 


2r 1 


(s) = r 4 — (\fl+e a^-1) + — (V'ea-L - In |\/ea^ + *s/l+e a ± 

O 3.-1 cL i \ ~ ~~ ~ 


r 2 2b n , 

+ r " h — - e 5 / 2 + 

o 2 5 € 


(4.27) 


and 


2x' 


x (s ) = x 4 — (Jl+ea^ - 1) ± — (\/ ea-i_ sfl+ea^ - In | \l ea~^_ + \ll+ea ± 

a i a i \ 

(r' + x' *>/ eaj ,(r"+b 1 e 1//2 )e \T l+ea x 




(x' ± r * */eab (x' ± r ' s/ea^) 


+ ...) (r"+b 1 e 1 / 2 )ede 


+ . . . . 


(4.28) 
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For the determination of the two constants a ± and b l5 we employ the points 
v = j 5 k, the two points immediately downstream of the separation point. 
With the abbreviation A = rl - - r'^e (l-(€j/e k ) 1 / 2 ), 


f k-T ’ 


= - <w 3/ii — L - 

J N/l+e k a/ 


Si “ 


1 +£ j a i 




\ 

i a V J " (£,/e k ) 3 ^ 2 ) j 

V/l+e a x J ^l+e k a i ' ^ 


Since a x appears also on the right side, we solve the expression for a x 
with the regula falsi . The constant b x is now^iven as 


bn = 


7372 

e k 


, I _ 


1 + x' VeTai 




The curvature of the free streamline can now approximately be written as 


+ x'n/ ea, a 1 (r ' + x' \/ ea x ) 

o 2 Z o 0 

2 \fe" (x^ ± r^ \l ea x ) 2(l+ea x )(x^ ± r^ «>/ ea x ) 
(r"-3b x e l/2 /2) v/l+ea x 

0 .... (4.29) 

(*o ± r^ 4e&l) 


(s \ = r"(s l = 
* (s) x'(s) 
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If the flow separates smoothly from a carved obstacle, the constant a x 
disappears and the singular term drops out of the expression. 

r" 

lim K(s) = . 

€-> 0 X 0 
a l-> 0 


This is the curvature of the contour at the separation point. 

For this case the derivative of the radius with respect to 
the arc length is, with one additional term, 

r 1 (s) = r^ + r^e + h-^ 3 / 2 4* r^ f € 2 /2. 

The second derivative is continuous since the singular term vanished. 
The third derivative, however. 


r m (s) 


2 Ll 


^ + 


can be singular, depending upon b x becoming zero or not. The derivative 
of the x-coordinate becomes now 


x 1 (s) = x f 


o o 


r"e 2 


€ + 


2x* 

o 


b r* 

+ -A 2 - e 3 / 2 + 


> r"e 5/2 
1 o 


X 


b ^€ 2 

2x f 

o 


where -x" = r^r^/x^ and K o = r o^ x o 9 To determine the constant b X9 we 
employ the point y =k , downstream of the separation point and retain 
only the term with b x : 


bl ■ [r k ■ <r o + r o 6 k )I/e k /2 - 


The radius of the free streamline becomes now for the smooth separation 


r(s) = r +r l e + r"^- + hi e5/2 > 
o o o 2 5 x 
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and the x- coordinate is 


x(s) = X_ + x 1 € 
x ' O O 


o o 

2x* 

o 


2b r 


>2 _ 


X o 


5x* 

o 


: 5/2 


r"e 3 

o 

6x f *** • 

o 


The function (4.24) is chosen more or less intuitively; only for the 
Helmholtz flow around a disk has (4.24) a justification to represent 
the free streamline in the vicinity of the separation point. However, 
equation (4.24) was fitted to a great number of calculated cases, all 
of which the curve fit was very good and represented the geometric 
functions of in the immediate neighborhood excellently for the 
abrupt as well as for the smooth separation. 

For the determination of certain velocity components, for 
instance, the tangential velocity, the series expansion of the Kernel 
function needs the derivative of the source strength, dq f (s)/ds. In 
general, the source strength varies with the derivative of the contour 
as indicated by equation (4.1). The derivative of the source strength 
behaves, consequently, like r n (s). We assume, therefore. 


q(s) = q o + Cx-fe + q^e + C 2 e 3 ^ 2 . 


The points v = j 5 k downstream of sg are again used to determine the two 
constants C x and C 2 . We designate the value of the source strength at 
s B as q Q . We obtain this value by a simple extrapolation of the source 
strength at points upstream of s^. It is 


and 


Ci = 


e x!i 




V L 


q, - (q +q* e, ) 

k o o k' 


-3/2 


(q +q* e . ) 
o ^o Y 


€ 3 


Co “ 


= q i ~ (q o +q o e i ) _ 


-3/2 


e . 
J 
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The derivative of the source strength with respect to the arc length is 
finally 


q ' ( S ) = i- + q' + C 2 f JT + ... . 

2 ve ° z 


Again, the constant C x disappears if the separation of the free stream- 
line is smooth. 


(4) Certain Changes of the General Procedure for the 

Application of Other Models 

In the general discussion of the theory, we did not mention 
the changes for other than the Helmholtz model. As we already mentioned, 
this model has infinite streamlines r&. The integration, however, is 
terminated in s E . This point must be placed sufficiently downstream of 
s E in order to render the errors small. A special error investigation 
approach is given in appendix D which determines, when applied, whether 
the termination point s E should be placed still further downstream. 

Levinson [17] determined the asymptotic form of the free stream- 
line. With the assumption of suitable regularity conditions, he showed 
that the free streamline had the asymptotic form of 


r 


c* 4T 

(In x) 1 / 



1 In (In x) 
8 In x 


(1/ In x) . 


The constant C* is a form factor and can be linked to the drag of the 
body by the relation 


D = f PC* 4 Uoo . 


Sufficiently downstream of Sg, the slope of the free streamline is small 
and we can assume, according to the slender body theory, that the source 
strength is proportional to the slope of the contour 
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1 


/• \ dr 
q o (x) 


2 ^ (In x) l/4 


1 - j (1/ln x) 1 / 2 + 


We observe that the source strength decreases in proportion to C*A /x~. 

The asymptotic form for q Q (x) was fitted to q Q (x) of table 2, and 
an excellent agreement was obtained. 

According to equation (D. 18) of appendix D, the induced 
normal velocity at s B of that part of the infinite wake which is 
neglected is proportional to the source strength at sg, to the inverse 
of x E , and to the inverse of the exponent (k+1), which is a measure for 
the decay of the source strength q with increasing s. Since in the 
Helmholtz model the source strength attenuates only moderately with x, 
the end point s E or x E has to be taken rather large in order to keep 
AV^(s b ) small. 

The Riabouchinsky model requires rather extensive changes in 
the method of solution. The symmetry of the model is reflected in the 
symmetry of the matrix elements: a v ^ - a ( m - y) ( m - > where m = N + 1 . 

Due to the model symmetry the calculation effort is reduced to half. 

The right side of the system of equations is an odd function (b v =-b( m _ v \) 
with respect to s R for 2r{, and (cp Q ) v , it is an even function (b v =b( m _ v \) 
for the constant C. From this behavior of the system, we can deduce that 
(q 0 T v is an odd function and (A^ 0 ) v is an even function with respect to 
sr. We have now the possibility to combine certain elements of the 
matrix. For the numerical solution, we use only the N/ 2 points on the 
forward position of the body and streamline up to the point of symmetry 
s R . 

We write the matrix as 


. (4.30) 


A 1 1 A 1 2 •* A IN/2 


q + Aq 

H oi H oi 


2v x + 0 

A. A. .. A , 


q . +Aq . 


2r . + 0 

11 12 lN/2 


Ol Ol 


1 

A. A. .. A . 


q . + Aq . 


cp + 1 

j 1 J 2 JN/2 

• • • • 


°J oj 


OJ 

V 21 \/2 2 •• ^/aN/a 


q oN/ 2 +Aq oN /2 


^ON /2 + 1 
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The elements are now formed by the difference or sum of two a v ^' s . 
the region 1 ^ v § i, the elements are formed by the expressions 


A = a + a 
vv vv v(m- v) 


1 + 


2icr 



8r 


-i 

"v 


* S v” S E ' 


+ s In 1 - — 
V 1 s 

V 



[ 1=1 



8r 


v 


s -s 


V 



r v K v 



r 

^ m- v 


*/2 


[r x ! -r f (x -x )]G (k 2 . N )-r x'G^k 2 *) 
1 v v v v v tn-y _ o v v(m-v) m-v_v. x v(m-v) 

[(x -X ) 2 + (r +r ) 2 ] 3 / 2 
v V m-v v m- ir 


L (m-v) 


m- v 


and 


r( r ,x'-r' (x,-x ))G rt (k2 )-r | r'G 1 (kf ;ii ) 


A = a + a / ~ 

v\i V[± v (m- |jl) jc 


V v v v v ii o yii . |i v 


via 


« vV ! + < r vV 2I3/2 


r A 
M- H 


- [r y x y~ r y (x y~ X m- u ) 3 G o (k y fa- uV > (m- u.)K G ^K(m ru ) } 

+ < r v +r (m-n )> 2]3/2 


r A 
m-v m-v 


For points on the free streamline (j ^ v ^ N/2), we obtain 


A = a + a 
vv vv 


= U ± 

v(m-y) it [2 


_ N , 



8r 


In 


i-l 


s -s 
H v 1 


8r 


s„ In 


s “ s 
1 y E 1 


+ s In 1 - — + 


r (m- y ) F o ^ k y (m- y) ^ ^m-y 

[(x -x ) 2 + (r +r ) 2 ]3/ 2 
v m-v v m-v 


For 
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and 


A 

vji 


i + a 
V[i v (m- \x) 


A r F (k 2 ) 

M- [X o VU 


(x -x ) 2 + (r +r ) 2 


+ 


A (m- u) r (m- a) F o y (m- u) ^ 


<7 (x -x ) 2 + (r +r . N ) 
v m-(j. ' v 



If we solve the system for the function (q o ) v > we apply the 
upper sign. The lower sign is used, however, when (Al 0 ) v is the solu- 
tion to the system. 

The procedure to solve the system (4.30) is the same as it is 
with the Helmholtz model. The position of the plane of symmetry is 
found with the condition r 1 (s) = 0 on r*. We call this point s r. 

The solution of the system (3.2) provides us with the perturba- 
tion potential (cp*) v and the tangential velocity (V$) v on the contour 
Ti and P %. As a next step, we determine with (4.8) the expression (cp Q ) v . 
The perturbation potential of the symmetric body-cavity configuration is 
an odd function with respect to Sr. However, it can only be an odd 
function if the integral expression in (4.8) disappears. In general, 
this will not be the case; therefore, we have to subtract from the func- 
tion (cp o y v the value of this function at Sr in order to make (cp Q ) v an 
odd function 


cp o (s) 


?>) + 


(V, 


TB 


V*(s)) ds 


c Po (S R ) ' 


'B 


In the past we could freely change the constant C. For the solution of 
this system, we again set C = 1. With equation (4.16), we determine the 
constant C. If the constant C does not vanish, the general solution 
(q Q ) v is not an odd function with respect to s R , since it still contains 
the even contribution C(£q Q ) v . The normal velocity (Vn) v is calculated 
along the entire length of p v , and with it and equations (4.17) and 
(4.18) an improved line 1^* is determined. Again we apply the condi- 
tion r 1 (s) = 0, and a new point of symmetry sr is obtained. The cycle 
can thus start anew. 
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Few changes are required for the dissipation model . For 
cavity pressures Cpg < 0, the absolute value of the curvature has a 
minimum downstream of sg; thereafter it increases again. The curvature 
increase of the free streamline in the theoretical flow will not be 
realized in the physical flow because of turbulent dissipation. The 
tangential velocity Vtb along the free streamline is constant only for 
a short distance behind the base of the body. After that it decreases 
in order to obtain the frees tream value at infinity; V-p = 1 . Therefore, 
the free streamline of the dissipation model is replaced by a straight 
line T 3 . The contact between r* and r 3 is of first order. r 3 repre- 
sents the trace of a stream tube of constant diameter. The point of 
contact is sp, where r* (s) = 0. The contact point s D is to be deter- 
mined for each iteration step. It represents also the end and starting 
points of two characteristic sections of the body. Therefore, the fixed 
points in these two sections separated by s^ must be determined for 
each step 


The stream tube r 3 is considered to be a solid surface. Con- 
sequently, the boundary condition (1.21) must be satisfied on r 3 . The 
integral equation (4.1) or (4.9) which is to be used for the region 
and which satisfies (1.21) is now also applied to the region r 3 . Thus, 
the linear algebraic equation has the following matrix form: 


a ll 

a 12 

• • 

a iN 


q + A} 

01 Ol 

• 

• 


• 


• 

a . 
ii 

a 

i 2 

- - 

a iN 


q . + Zfl . 

01 01 

a . 

a . 

J 2 

• • 

a jN 


q . + Aq . 

Oj ^Oj 

• 

• 

• • 

• 


• 

a ki 

a k 2 

. - - 

a kN 


q ok + Aq ok 

a li 

a i 2 

• • 

a iN 


%r +Aq oi 

• 

• 

9 • 

• 


* 

a Ni 

a N 2 

• • 

a NN 


q oN +Aq oN 








_ _ 


2r , 1 + 0 


2r * + 0 
1 


CD . + 1 


^ok 


0 


0 


+ 1 


+ 0 


+ 0 


(4.31) 


c 0 



The elements of the matrix (4.31) for the regions 1 g v < i 
and j ^ v ^ k are exactly the same as those of the system (4.11) for 
the regions 1 g v ^ i and j ^ v ^ N. The elements of the additional 
region r 3 (£ ^ 0 ^ N) of (4.31) are given by the relation 


a 


vv 



r vV E 



and 


a 

vp. 


r (r x 1 -r 1 [x -x ])G (k 2 )-r x , G 1 (k 2 ) 
J± V V V v u o vu u v 1V v\x 


«/2 


[(x -x ) 2 + (r +r ) 2 ] 3 / 2 
V V \x V V \±' J 


A . 
M- 


The right-hand sides of (4.31) for the same region are 


b = 2r f = 0 
v v 


if we solve for (q Q ) v , an d 
b y = C = 0 


if we solve the system (4.31) for (£q Q ) v . The source strength q Q (s) is 
continuous at the point of contact s^; its derivative, however, is dis- 
continuous and we have to treat q Q (s) in the neighborhood of sp separa- 
tely. The integral of (4.8) is evaluated only between Sg ^ s ^ s^, 
since (cp Q ) v is used only along r%. 

The position of Sj) is a function of the cavitation coeffic- 
ient Cpb* With decreasing Cpg, the point of contact approaches the 
separation point Sg. For the numerical calculation, we terminate the 
stream tube r 3 a certain distance sp downstream of the separation point. 
According to the error estimation given in the appendix, it is sufficient 
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to let F 3 extend approximately 10 radii downstream, since q 0 (s E ) attenu- 
ates rather quickly on r 3 and the exponent k is rather large. 

For the investigation of the rearward separation point from 
a sphere, the free streamline has to be altered in the vicinity of 
the x-axis. We understand by the rearward separation points those 
points sg for which the free streamline forms a cusp or intersects the 
x-axis. For these cases, the curvature of the free streamline is posi- 
tive, and therefore concave to the outer flow field. Due to the vicinity 
of the x-axis, the nature of the series expansion (1/r » 1), and the 
inadequate treatment of the rearward stagnation point, small errors in 
q Q (s) appear in that region. The normal velocities computed with these 
q Q (s) tended to make T% convex; whereas, the real solution should remain 
concave in order to form a cusp. Therefore, r% was replaced in the 
region of the rearward stagnation point. Whenever the curvature of the 
free streamline became lower than a certain value k q , the curvature 
K( s) was replaced from this point on with an exponential function 


K(s) = K 0 exp(-a(s-s L )). 


The factor a was determined thus: that the next derivative of the curva- 

ture K 1 (sg) was continuous. The curvature k q = k( Sg) can be determined 
with an expression derived from (C.16) by requiring that the tangential 
velocity at sg be equal to the specified velocity Vgg. 

We write 


V t k = c>V T /c>v = (l-V T )/^h 


with an unspecified Z^h. Next, we form the ratio k q / K /( ( s^) and set 
Z^h’Vztfi = 1. Thus, we obtain an approximate equation to determine k q 


K o = 


K < S L > 


v t( s l> C 1 - 


V TB^ 


V TB^ 


V£(s l )) 


TB’ 


vv 


* 1 


where K* (sg) is the curvature of r% at Sg. On the remainder of the line, 
which we will designate r 3 in accordance with the dissipation model, the 
pressure cannot be specified any more. Therefore, the rearward stagna- 
tion point location is in error, and also the curvature at that point is 
wrong. Another method to find the cusped cavity is mentioned in a later 
chapter. 
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The system of equations for determining the source strength 
is the same as given by (4.31) with the exceptions of the right side 
for region £ ^ v i N, which becomes now 


b v - 2r ^ 


if we solve the system for (q 0 ) v , and 


b v = e = ° 


if we solve for (£q 0 ) v . Since the derivative of the curvature is con- 
tinuous in s^, the derivative of the source strength dq Q (s)/ds is con- 
tinuous and smooth in the neighborhood of that point. Therefore, no 
special treatment of the region adjoining s-^ is given. 

For a certain (s^; Cpg) combination with Cpp > 0 and C D < °> 
the free streamline does not intersect the axis of symmetry. The 
curvature of the free streamline is still positive and concave; how- 
ever, it has a minimum downstream of sp and after that, it increases 
rather rapidly. The function r'(s), originally negative, changes its 
sign in the vicinity of the curvature minimum and r 2 diverges. Whenever 
this happens, we determine the point of contact between an d F^ with 
the condition r* (s) = 0 and apply the concept of the dissipation model. 
With Cpp decreasing and Cj) approaching zero, the stream cylinder r 3 will 
diminish until finally Cj) = 0 and r 2 is forming a cusped cavity. Thus, 
an extrapolation of the above-mentioned dissipation model will lead to 
the cusped cavity. 

In the general discussion, it is implied that the parameter sg, 
the separation point, is given and it is kept constant during the itera- 
tion. However, if one wants to determine the point of smooth separation 
from round bodies, one must apply an iterative procedure. For a number 
of tentatively chosen separation points, one conducts the complete iter- 
ation procedure as described and determines the constant |\fa^| of equa- 
tion (4.24). If two points sg are encountered between which the 
expression ± \Ta^ changes sign, one can easily determine the point for 
which a-L vanishes or becomes sufficiently small. This is the point of 
smooth separation Sg (smooth). By plotting the curvature of F^ versus 
(s-Sg) logarithmically, one can decide if Sg can be improved. At the 
same time, one can check (4.29) for its validity. 


53 



The Cavity of the Lifting Body 


The cavity of a body placed at an angle a toward the flow will shift 
to an asymmetric position. The body experiences a lift proportional to 
the angle of attack a . The circumferential velocity Wx(s,q:) on the body 
surface will also exist on the wake or cavity surface, thus creating 
locally over- or under-pressures which in turn will reshape the cavity 
and will destroy the rotational symmetry. The local pressures will shift 
the wake surface to such a position that the condition of pressure con- 
tinuity across the wake surface is satisfied. We assume, again, negli- 
gible velocities within the wake. Consequently, a necessary conclusion 
is the vanishing of the local normal force along the wake or cavity. 


C n (s,a) = 0. 


(5.1) 


Our main concern in this report is not so much the proper conditions of 
the flow along and around the wake stream surface and its exact position 
in space, but the influence of the wake on the lift distribution of the 
forebody. An exact treatment of the wake especially of the distant 
wake is not necessary since velocity and thus pressure perturbations 
attenuate rather quickly. We employ therefore a simplified model. Its 
schematic is given in Figure lc . Since we treat again only bodies with 
rotational symmetry, we must assume that the flow separates along a cir- 
cumferential line which lies entirely in one plane normal to the body 
axis. This exemption excludes round bodies with an interior separation 
line . 


Condition (5.1) requires a certain tangential velocity. We can 
calculate it and call it With this velocity we can define a 

potential 0 lB (s) which is known along except for a constant C. 

To describe the general theory we employ again the Helmholtz model 
(or infinite streamline model). We assume for reason of simplicity that 
the wake pressure coefficient Cpg is independent of the angle of attack 
qj. This is a reasonable assumption for small angles of attack. We 
therefore use the same Cpg of the axisymmetr ic flow case. For the 
cavity subjected to under-pressure we will use the dissipation model 
exclusively. Special changes of the general theory will be discussed 
at a later time. 
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(1) The Potential and the Velocity Components 

We place now the forebody at an angle of attack a toward the 
parallel flow field. The wake will shift to a new equilibrium position 
where condition (5.1) is satisfied. The exact potential of this body 
will be developed now. Assumptions (1) through (3) of chapter 4 are 
also valid for this case. 

The total potential of the body is in body-fixed coordinates, 
x°, y°, z°, with the assumption that the x°,z°-plane contains the velocity 
vector given as 

01 (x°, y°, z°) = z° sin a + cpxCx 0 , y°, z°) . (5.2) 


The axial component of the undisturbed flow is sin a, the normal com- 
ponent is Vcq = Uoo sin O'. Only this component will be considered here. 
The first term on the right-hand side of (5.2) represents the potential 
of this field. The function cp 1 (x°, y°, z°) is the perturbation poten- 
tial of the body with bent cavity. 


<M*°. y°, .») - £ ff , 

^ n / ( x °- i °) 2 + ( y °- T ) 0 ) 2 + ( z °-^°) 2 


(5.3) 


where dS designates a surface element. 

The potential (5.2) is for our purpose still too general. We 
therefore want to make it more accessible to calculation by using certain 
characteristics of the cavity at small angles of attack. Since we deal 
usually with the pressure distribution of the forebody an exact descrip- 
tion of the wake is not necessary. The only M lif t-producing M component 
along the wake is sin (o'(s)). This component will shift the cavity 
centerline, which can now be described by the function z 1 (s) of Fig- 
ure lc . Our goal is to find the position of the centerline of the cavity 
or wake which will be attained due to condition (5.1). 

To formulate the problem, we make the following simplifying 
assumptions : 

(a) The cavity cross section remains rotationally symmetric. 
Its diameter is fixed and obtained from the axisymmetric 
flow case. 
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(b) The planes normal to the body axis are shifted parallel 
to itself. The new trace of the bent axis in this plane 
is situated along a line z x (s). 

(c) The local normal force vanishes along the wake or cavity. 
This requires a special function z^s). The velocity 
V^b known as the cavity streamline; so in the total 
potential jZTjjB (s ) except for a constant. 

(d) We neglect locally the condition of pressure continuity 
across the stream surface. The condition is satisfied 
only in the mean according to assumption (c) . 

(e) Only terms linear in a are considered, higher terms are 
neglected. The contour of the cavity or the wake is 
assumed to be small and negligible, 8(s) « 1.0. 

(f) Only those body forms shall be considered here which 
have the line of flow separation in one plane normal to 
to the body axis. This excludes bodies which permit 

an interior separation. 

(g) Since the exact angle of attack distribution a(s) along . 
the free streamline is not known a priori an approxi- 
mate distribution QT f (s) must be assumed. The resulting 
cavity boundary is designated as F%. All calculations 
are performed approximately with cn*(s) . 

(h) Normal velocities (Vni(s,QI) £ 0) will result on r* if 
a* (s) and a(s) are not identical. The magnitude of the 
normal velocity indicates how and how much the angle of 
attack distribution must be changed. 

The next step is the development of the potential in cylindrical coordin- 
ates. We refer to Figure lc where the coordinates of a fixed point po 
on the surface of the body according to the assumptions are given in 
body-fixed coordinates as : 

x° « x; y° « r sin to; z° ^ z ± (s) + r cos go; s° = s. (5-4) 

A simple integration from the separation point sg to a field 
point s yields the function by which the wake is shifted 



Since the arc length is approximately the same as in the axisymmetr ic 
case we keep s as the independent variable. For the summing point Q° 
on the surface of the wake one obtains similar to (5.4) 

6° « T)° « p sin u) 1 , £ 0 » z x (a) 4- p cos co 1 ; a° « a. (5.6) 


The potential of the undisturbed flow field, represented by 
the first term on the right side of equation (5.2), can now be written 
as r(s) cos u) sin O'. The denominator of the integral expression (5.3) 
contains the distance between point P° and Q°: 


R = n/ ( x °-^) 2 + (y°-T] 0 ) 2 + (z°-^°) 2 . 


(5.7) 


If we introduce (5.4) and (5.6) into the square root expression for R, 
we obtain finally for 


1 = 1 

R (x- £) 2 -fr 2 4-p 2 - 2r p cos (w-co f )+2 (z x (s) -z !(a)> (r cos co-p cos to 1 )*f0(5 2 ) 

(5.8) 

The relation between the angle 6(s) which we call camber and the angle 
of attack a(s) can be seen in Figure lc : 


a - a(s) + 5 (s) . 


(5.9) 


Consequently, the function z x (s) in (5.5) becomes 


s 



a J 


ds . 


(5.10) 


For small values of the angle of attack a, one can develop the recipro- 
cal distance R into a series which we will truncate after the second 
term: 
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1 


1 

R 


n/ ( x- |) 2 +r 2 +p 2 -2rP cos (to* -w) 


2q(zi(s)-zi(g)) (r cos pj-p coscj' )/g 
[ (x- |) 2 +r 2 +p 2 -2rp cos (w 1 -oo) ] 3 / 2 


(5.11) 


For the source distribution of the body and the cavity for a small angle 
of attack, we choose the symbol (i(s,a) = q 1 (s)a(s)/a. After we insert 
|i(s,a) and equation (5.11) into the perturbation potential (5.3), we 
obtain the approximate potential of the body with a cambered cavity in 
normal flo^. 

S E 2rt 

<Pi(s,a,w) « - ^ J J 

0 0 


^ g ^ _ fi . cos_. M ' _ . 0(a2 A daj , da< 

(v-n 2 4-r 2 4-A 2 -9rn / l \ ' 


The comparison of this expression with equation (1.1) reveals 
that the first term under the untegral is equal to the perturbation poten- 
tial of the body with a straight cavity. The next higher terms are at 
least proportional to a 2 . If we apply partial differentiation and obtain 
the velocity terms depending on a 2 , we can prove that these terms do 
not contribute to the total lift, since the exponent of cos co is even, 
and the integral of the expression for the normal force vanishes. We 
consider, therefore, from now on, only the perturbation potential 


9 1 (s,q, w) = - ^ 


E 2jt 

a P P 4(q,cn) p cos q/dc o 1 da 


0 0 


v (x- £) 2 +r 2 +p 2 -2r p cos (co 1 -gj) 


which is the potential of the axisymmetr ic forebody-cavity combination 
with a variable angle of attack distribution along the wake line r|. 
The total potential is now given as 


0i( s >Q!,gj ) = r (s ) cos co sin a + cpi(s ,a, go) . (5.12) 


On the surface of the fixed forebody, the angle of attack is 
constant and we set sin a = 1 . Along the angle of attack a(s) is 
determined by (5.9). The angle a is still a constant; however, it con- 
sists now of two variable contributions. 
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Partial differentiation with respect to the arc length s 
yields the tangential velocity for meridian plane co = 0: 


V Tl (s,a) = (sin a + v x (s,a)) r 1 (s) + u 1 (s,a) x f (s). (5.13) 

The normal velocity for the same plane is given by 

V Nl (s,oO = (sin a + v 1 (s,a)) x* (s) - u 1 (s,a) r' (s) (5.14) 


and obtained by partial differentiation with respect to v. 

These velocity components as written above are valid and 
exact only along r l5 the contour of the forebody. For obtaining the 
velocity components on the cambered cavity we have to insert (5.9) 
into (5.13) and (5.14), linearize according to assumption (4) and 
the resulting approximation velocities are: 


V Tl (s,a) = (a(s) +v 1 (s,a)) r' (s) + u 1 (s,Q') x f (s) (5.15) 

and 

V Nl (s,a) = (QJ(s) + v 1 (s,a)) x’(s) - U 1 (s,a) r f (s). (5.16) 


The terms Q'(s)r* (s) and a(s)x* (s) are in this form the tangen- 
tial and normal velocity components of the normal velocity u ro sin (a(s)) 
of the undisturbed flow field respectively. The perturbation velocities 
u i(s,Q!) and v 1 (s,-a) are those of the straight cavity with variable angle 
of attack distribution. 

The linearized problem consists now in finding the angle of 
attack distribution a(s) along the straight and rotational symmetric 
wake or cavity which produces a vanishing normal force C n (s,oO = 0 
along the line r%. 


(2) The System of Integral Equations 


An integral equation with which we are able to determine the 
source strength on is given by equation (5.14) with V]S[ 1 (s ,a) = 0. 

We introduce the perturbation velocities, linearize, and obtain 
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(5.17) 


n(s,a) = - 202 c' (s) - -J 1 (i(a,a) K q (s,a) da. 


The Kernel is represented by equation (1.28), where n = 1. 

For the cavity boundary r%, we use the total potential 0j^(s y a) 
expressed in terms of the perturbation potential with unknown source 
strength. With assumption (c) the potential 0jj&(s,O0 is known. With 
equation (5.12) for the meridian plane co = 0 and a replaced by of (s), the 
assumed angle of attack distribution, the integral equation for the 
region r* can he written then as 


0ib( s >°O - r' (s) a*(s) 


-hf 


M.(cr,a) 


K 9 (s,a) 


ds . 


(5.18) 


where the Kernel function is in accordance with equation (1.9) and n = 1 
given as 


V s 5 cr) 


-2pFi(k 2 ) _ 

n/(x -^) 2 + (r+p ) 2 


In both integral equations, (5.17) and (5.18), all functions except the 
source strength [a(s,a) are known or assumed to be known. The Kernel 
functions depend only on the geometry of the forebody and the cavity. 

We employ according to assumption (a) the geometry of the axial sym- 
metric flow case. Thus, Kq(s,a) and Ky(s,o) do not change any more, 
and we need to determine them only once at the beginning of the itera- 
tion. The angle of attack along r x is constant and a = 1. Along 
however, a(s) must- be assumed first. 

On the left side of equation (5.18) the potential # 33 ( 8 , a) has 
to be evaluated. The next step toward a solution is therefore to relate 
the potential to the normal force coefficient. We start out from the 
definition of the local normal force coefficient: 
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(5.19) 


C n ( s ,CC) 


-r (s) 



(s,a, to) 


O 


COS OJ d GO. 


The pressure coefficient again can be expressed by the tangential veloc- 
ity. In conjunction with the velocities of the axial flow case, we 
obtain the total tangential velocity as 


V^ G ( s .«> w ) = (V T (s) cos a + V Tl (s,a) sin a cos to) 2 
+ v| 2 (s ,a) Sin 2 a s in 2 to. 


(5.20) 


where the velocity component in circumferential direction is given as 


v T2 (s ,a) = (-sin a(s) + w 1 (s,a)) . 


Since in the equations the velocities V«p (s,o;) and Vq* 2 (s,oj) are con- 
nected with the angle of attack, these functions are therefore based 
on sin OC = 1. The perturbation velocity w-lCSjQO is given that value 
which it will attain in the meridian plane jt/2. For a positive angle of 
attack, this value will, in general, be negative. With the retention 
of all terms, the pressure coefficient becomes 


Cp (s ,0!, to) = Cp (s ,Q!=0) 


\/ 1 - C p (s ,o:=0) 


V t 1 ( S ,Q; ) S:i ' n COS w 


+ (1-Cp (s ,Q!=0) - V 2 i (s,a) cos 2 co - V T2 (s,q:) sin 2 co) sin 2 o;. 

(5.21) 

Insertion of this expression into (5.19) produces no contribution of the 
first term, and third (a 2 ) term, only the term linear in a contributes 
to the normal force; hence 


c n (s,Q!) = 2a: it r (s) 4l - Cp(s,a=0) v Ti (s,o:). 


(5.22) 
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The integration of V 1 j, i (s,q:) along s yields 


0 i (s, a) 



c (s,a) 

° — ds + 0i(s a). 

r(s) a v l-Cp(s,ce=0) 


(5.23) 


The normal force shall vanish along the exact free streamline r 2 , hence 
the potential becomes 


0±B< 8,a ) = 


(5.24) 


The cavity boundary is thus an equipotential line, which means that the 
tangential velocity vanishes on Fg. We set Vrp 1 (s,o;) = 0, also V^ i (s,cc) = 0, 
since F^ i- s a streamline and we arrive at a condition that 


u i( a >oD = Ot(s) + v 1 (s,a) = 0. 


(5.25) 


The constant ^(s^a) in equations (5.23) and (5.24) is the value of the 
total potential at the separation point s B . 

The left side of the integral equation (5.18) contains the 
expression (j2fiB (s ,Q!) - r(s) d*(s)) which could be construed as the per- 
turbation potential of the cavity in normal flow, if a ,v (s) would be 
identical witha(s). However, we are using for the numerical calcula- 
tion the assumed angle of attack distribution and therefore 0-j} (s, a) - 
r(s) a*(s) is only approximately equal to cp 1 B(s ,a) . For a convergent 
procedure F% F 2 and Q? v (s) ->Q:(s) and, consequently, we define 


cpiB( s ,a) = ^33(3, a) - r (s) ct (s) = cp 1 (s,a) + c. 
The barred term represents all known functions 


cpi(s,a) = -r(s) a*(s), 


(5.26) 


(5.27) 
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and the constant becomes 


c = 0iB( s B ’ a )* 


Similar to the axial flow case, we can represent the perturbation poten- 
tial of (5.26) with the perturbation potential cp*(s,a) of the assumed 
angle of attack distribution plus some corrective terms. After combin- 
ing 


cpi(s>a) 


0t(s B ,O!) 


- r (s) cf(s) + 


s 

I 


v* x (8,a) 


ds 


B 


and 


^3J5 (S B ,a) " 0l( s B ,a) = C > 


we add and subtract in equation (5.26) the terms 


0i(s B ,a) 


and 


/ 


v^Cs.a) ds 


to obtain finally, with equation (5.23), 


* 1 P c n ( s >a) ds 

cp 1B (s,0!) = cpiO.a) - 2^ / - 1 ■" + c - 

«T r (s ) n/1—C_ (s , qj=0) 


(5.28) 


The constant C depends again only on the potential difference at the 
separation point s^. The perturbation potential cp*(s,o:) and the normal 
force coefficient C§(s,qi) of the approximate free streamline T% are 
obtained with the solution p,(s,o:) of the integral equation (3.1) with 
n = 1 and q 1 (s,a) = -x(s)o£*(s). Since (3.1) was obtained with the 
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condition (s = 0 along the entire boundary r l5 r% and r 3 , F% is 
consequently for this case a streamline. We combine again all known 
functions and obtain 


cpi(s,a) = <p*(s,a) - 


C^(s,a) ds 


r(s) n/1-C p ( s,Ctf=0) 


(5.29) 


If the assumed angle of attack distribution Of(s) is approaching the 
exact one a(s), the integral* in equation (5.29) and the constant C in 
equation (5.28) will disappear. For the evaluation of the integral in 
(5.28) the integrated Lagrange interpolation polynomial with n = 4 
points, two on each side of the intervals, was used. 

The representation of the potential cpxBCSjO:) with the assumed 
perturbation potential cp*(s,a) and additional correction terms had 
advantages for the numerical solution of the integral equations. The 
discussion of these advantages follows the same line as that of chap- 
ter 4 of the zero angle of attack case. A repetition is therefore not 
intended here. 

Next, we continue with the discussion of the procedure to 
obtain a solution of the integral equations for the lifting body. 

Equation (5.28) is inserted into the left-hand side of the 
integral equation (5.18). For determining the constant C, one other 
equation must be furnished; temporarily we might set C = 1 for the 
numerical calculation. 

The system of integral equations (5.17) and (5.18) is converted 
to a set of N linear algebraic equations which are satisfied in N points 
on the contours F i_ and r%. We introduce summation signs instead of the 
integral signs and obtain now for the integral equation (5.17) of the 
region r 1 with points 1 ^ v ^ i 



r K 
v v 


A 


+ 



8 r 


In 


S -£ 
1 V 


+ s In 1 



v 



r K s 
v v 



(equation (5.30) continued on next page) 
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+ 


- y’ 

2n L 


i=l 


[r x'-r'(x -x )]G (k 2 )-r x' (G (k 2 )+G (k 2 ))/2 
2 v v v v u. 1 v|i |1 v v o v vn' g v V|i A 

^ [(k -x ) 2 + (r +r ) 2 ] 3 / 2 V 

V |J. V |i 


(5.30) 


Along Ti, the angle of attack distribution is constant, and oc v - 1 . The 
primed summation sign indicates that the summation at the point jj. = v 
is left out. 

For the region r% with points j ^ y ^ N, the following equa- 
tion determines the source distribution jj, y . The constant C, as men- 
tioned above, is set equal to unity. 


(cpl) v +1 = 2^ 


S' 

LL=1 


8r \ r 8r 

In -| ¥—r - 2 ) A - is In -i 

s -s_. | J |i I E I s. -s 


(1 v 1 


V_ + s ln h . _E| 


v E 1 


- s. 


1 Y ^u 2r u Fl(k vu )A U 


2rt 


i [ (x -x ) 2 + (r +r ) 2 ] 1 / 2 
i=l v u. V 11 


(5.31) 


The two equations (5.30) and (5.31) are solved for the source 
strength according to a special procedure. In matrix notation, we 
obtain 


a ll 

a 12 

a . 

a . 

ii 

12 

a . 

a . 

Ji 

J 2 

a Nx 

a„ T 

N2 


a fN 


iN 


a jN 


NN 


(11 + All! 


H. + An. 


^ + 


% + ^ 


-2x! + 0 


-2x'. + 0 

L 


9 


lJ 


+ 1 


'PxN 


+ 1 


(5.32) 
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The elements are formed by the right sides of equations (5.30) and 
(5.31). In particular, we obtain: 


a vv “ 


1 + 


dr f ; C 1 - 


8r 


l%" s E l v 


ln 1 “ ~ ~ 2s^ ) - r K s 

i ~ 1 ^ ' v v E 


N , 


M- ! 


.=1 


8r 


x' In -I 

v V s -s 

'(IV 1 


- 3 ) r K At 
/ v vj pj_ 


for the region 1 g v S i; 


r A [r x'-r'(x -x )]Gi(k 2 )-r x' (G (k 2 )+G 2 (k 2 ))/2 
a = ^ <■* v v v v l-i/ 1V ml u_v \ ° vu 2V vu J _ 

V|J " /2 frv + r ,.) 2 ] 372 


V 


for the region 1 g v S i and v ^ p; 


N , 


w 2n 


8r 


In 


H=1 


I V s vl 


2 ) A - ( s In 

H V E is -s 


8r s_ 

— — r + s In j 1 1 - s 

v "E 1 v V E 


for the region j g v § N; and 


M- 


r F (k 2 ) 

n i v m_ 


71 [ (x — x ) s + (r +r ) 2 ]i/ 2 

L v v p/ v v p' J 


for the region j g v § N and v t p. 

We solve the system (5.32) twice in order to obtain the 
particular contributions. For the first solution, we employ the right 
sides 


b 


v 


-2x' 

v 


for the region 1 g v S i with CC y = l , and 
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b v ~ 9 iv 

for the region j ^ y ^ N, where cp lv is given by equation (5.29). The 
particular solution of the system yields the source strength jj, v . For 
obtaining the second part of the solution, the right sides of the system 
are altered. The elements a v ^ are the same as for the first part. 
Analogous to the zero angle of attack case, the right sides of the sys- 
tem become now 


for the region 1 ^ v ^ i and 

b = 1 
v 

for the region j ^ v ^ N. 

The particular solution of this system is A| x v , which is 
normalized to one. Since the geometry of the free streamline of 
the axisymmetric flow case does not change any more, the elements of 
the matrix which consist only of geometric functions will remain the 
same throughout the iteration. The same is true for the second right 
sides (0; 1). Therefore, the second system needs to be solved only 
once. The solution Api v can be stored and used for the iterations. 

After the particular solutions have been obtained, we com- 
bine them to form the general solution with an arbitrary constant C. 

H v = + c (5.33) 


However, before we discuss the determination of the constant C, we want 
to write the summation expressions for the tangential velocity. We use 
equation (2.10), replace the integral sign by a summation sign, and 
paying attention to equation (4.13) for the tangential velocity, obtain 


(^Ti^v r v a 


V V 




(equation (5.34) continued on next page) 
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2r 


s ' 2r 
v v 


8r. 


s -s 

1 V p 


N 


i v'« [x ; (x v- x u ) - r /; ]G i (k ^-v r v (G o( k v^ 2 ( k vu ))/ 2 

"2it Z, : A • 


n-l ^ [( x v"V 2 + ( r v + V 2]3/2 


(5.34) 


The general solution (5.33) is now inserted into expression (5.34). We 
set (Vpj) = (Vp x g) = 0, choose the special point for which this equation 
is to be satisfied at v = j, which is immediately downstream of s^ and 
solve for C: 


N , 


C = 


li'.A. - n . Fo - ) u K . A - a.r 1 . 
J 2 Zj P JP P J J 

n-1 

N 

) An K. A + All.Fo “ Au'.A. 

Zj P JP p 2 “ J 

u=l 


(5.35) 


with the special functions 



and 


K. = 2r 
JP P 


[x, i (x r x p) +r .i r j ]G ^ (k jp)- r P r j [ G o( k P^ (k P ]/ 2 

[ (x -x.) 2 + ( r +r.) 2 ] 3//2 
p J P J 
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Once the constant C is known, we insert the general solution into the 
expression for the normal velocity (2.9), which is now 


< V Hl>v 


a x 


v v 



1 




r K 
v v 


} 


A 

M- 


1 

2ot v 



8r 


IV S E I 


4- s In 
v 


1 




r 


K S 
V V 




2[i r 


[ r v x k r v^y- x ^J G i( k p- r /^ G o< k v l i) 4 <;g ( k v t i)J /2 

[(x -X ) 2 + (r +r )2]3/2 
V JJL V [X 


A . 


(5.36) 

The boundary condition V^(s,q:) = 0 is not satisfied on 
until the iteration cycle converges within a prescribed tolerance. 

Then [_i(s,ai) is also a solution of (3.1), and C n (s,cO vanishes on the 
cavity boundary. In analogy to the axisymmetric flow case, one can 
construct, with the aid of the velocity diagram of Figure lb, an 
approximate expression for the normal velocity on r%: 


v Nl (s,a) = v T (s)(e**(s,a) - e*(s,a)), (5.37) 

where 0 xx (s,a) is the new improved slope of the contour r^' • This 
slope can be approximated by the local angle of attack 


e** (s , a) = a - a** (s ) . 


(5.38) 


We insert this expression into (5.37) and obtain a new angle of attack 
distribution along the cavity or wake; 


a** (s ) ssq:*(s) 


v Ni (s,q:) 

V T (s) 


(5.39) 
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This distribution forms the basis for the next iteration step. The new 
angle of attack distribution CT*(s) is inserted into equation (3.1) which 
is then solved for p**(s,a). With this source strength, we determine a 
new normal force coefficient Cft* (s ,a) on r**. The perturbation potential 
(5.29) can now be improved. For this improvement we need to solve the 
system (5.32) for p(SjQ')* This procedure is repeated until condition 
(5.1) is sufficiently satisfied and the constant C is small enough. 

At the beginning of the discussion, we mentioned that the 
representation of the potential in terms of the perturbation potential 
of the assumed line Y% plus corrective terms received preference over 
any other approach. According to the calculated cases, we found out 
that the matrix is not as well conditioned as that of the zero angle of 
attack case. The jump of the function (0, 1) of the right side again 
caused fluctuations of the solution ( Ap x ) v in the neighborhood of the 
separation point s-g. These fluctuations entered the general solution 
with finite constant in all its severity. However, if the approach 
with vanishing constant C was used, these fluctuations vanished with C. 

For those cases for which the cavity is subjected to an under- 
pressure, we employ exclusively the dissipation model for the angle of 
attack case. In the axisymme tr ic flow case, we treat the boundary 
as a rigid stream tube with Vj^(s,a) = 0. For the angle of attack case, 
however, will shift and attain a position in space to satisfy condi- 
tion (5.1). Thus, it does not differ from fVp. Both r 2 and T 3 are 
therefore subjected to a varying angle of attack distribution a(s ) . 


(3) The Direct Solution 


of the Problems by Specifying V^B on r 


* 

2 


For 1 the case under consideration, the second approach is 
definitely possible by prescribing the tangential velocity V TlB on the 
cavity boundary as indicated in the previous chapter. For the region 
I* enclosing the boundary points j ^ v ^ N, we replace (5.31) by (5.34): 


( V TiR^ v 




In 



A 


P 



2r ) + 2r 
v / 


In 


8r -1 

l s v' s E l)}- 



2 m r 

p 


[< x 'v^v- x u )+r v r ; ]G i( k v U > - r tl r U G o( k ^- H; a ( k v| J >l / ? 

• [<* v -V 2 + <r v + V 213/2 . 


70 



The unknown function of the integro differential equation is 
the source intensity p v whose derivative is also sought. The geometry 
of the body-cavity configuration of the axisymmetr ic flow case is again 
employed here. The derivative of the source strength is obtained with 
a differentiated Lagrange interpolation formula [ 18 ]: 

p' = L' p + L'p + L' p . 
v- r v - 1 v*v v+i v+i 


The elements of the matrix for the region with the boundary points 
j ^ v ^ N are now given by 



and for v ^ p and p ^ v ± 1, the elements are 


a 

vp 


A r 
. IL-H . 

TC 


(x* {x - 

V V 


X u } ~ r v r P Gl(k vp ) " r p r v (G ° (k v^ )+Ga(k p )/2 

[(X v -X |a )2 + ( r v +r |i) 2 l 3 / 2 


The elements neighboring the principal diagonal contain, 
however, one additional term. It is the differentiated coefficient 
of the Lagrange interpolation formula 


a* , x = L 1 , + a , , v . 

v(v±l) V±1 v(v±l) 

This method was used for different numerical examples. The solution 
p(s,a) showed extremely severe fluctuations; thus, the normal velocity 
on r| could not be used to improve the angle of attack distribution. 

The procedure diverged, and it was therefore dropped. Even an increase 
in significant figures did not improve the situation. 
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The Forces on the Body 


To determine the total drag coefficient of the body, the component 
of the pressure coefficient in x-direction is integrated along the 
meridian angle go and along the body contour r x up to the separation 
point Sg. In addition, one has to add the base pressure coefficient 
which is considered to be constant over the entire base area. If we 
choose the base as the reference area, the expression for the total 
drag coefficient becomes 

S B 2rt 

C D CT m ~ J f c p ( s ) r ( s ) r ' ( s ) dojds " r 2 ( s g) 71 C p B’ ( 6 . 1 ) 

o o 


where cr m = R 2 jt is the dimensionless reference area. 

The drag of a body can also be determined with the cavity radius 
R" of the dissipation model. The radius R" is identical with the radius 
of the stream tube r^. We consider the body-wake configuration as a 
half body of radius R' r . We ask for the pressure drag of this half-body 
in potential flow. The drag is undetermined as long as the base pressure 
is not specified. We assume for this purpose that the half-body contour 
is interrupted by a slot in which the base pressure acts. Since the 
x-component of the integrated pressure along the contour disappears, 
only the suction force of the wake remains. We therefore obtain 


C 


P 


a. 2 



1 ). 


( 6 . 2 ) 


Additional forces and moments are the total normal force, the 
induced drag, and the moment coefficient. 

The force normal to the body axis, the normal force coefficient, 
is obtained by integrating the r-component of the pressure coefficient 
along the meridian and arc length direction; 

® g 2jr 

a C N°M = ~S r ^ x ' J G p( s ’ CJ » o: ) cos udu ds - (6.3) 

o o 
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The major contribution to the moment coefficient is obtained by 
multiplying the above equation (6.3) by the moment arm about a refer- 
ence point x(s): 


a " " 


r (s)x' (s)x(s) 


2 rt 

f C p (s, w,Ct:) cos co dco ds, 


(6.4) 


One other contribution to the moment coefficient is obtained by 
integrating the angle of attack dependent pressure coefficient about 
the body. In the angle of attack case, Cp is a function of the meridian 
angle go. The pressure coefficient therefore changes the sign with 
increasing go, and a force couple exists. 



r (s)r' (s) C n (s) ds, 


(6.5) 


where C n (s) is given by equation (5.22). 

Converting equation (5.22) to the pressure coefficient form, one 
term proportional to a 2 will appear which gives one contribution to the 
tangential force: 


1 


C T°M 


s 


B 



s )r 1 (s) 


o 


{2jt(l-Cp (ctf=0) ) - rt[(u 1 cos 0 + (1 + v-l) sin e) 2 


+ (-1 + Wx) 2 ] } ds . 


( 6 . 6 ) 


Considering all the contributions to the particular total force 
coefficients, one obtains for the induced drag of the body at angle 
of attack 


c D (aO = c N ^(sin a cos a) sin a + c T (a) cos a. 
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and for the lift coefficient 

C T (°0 = Ctj (sin q: cos a ) cos a - C m (o0 sin a, 
h T 

where the lift coefficient per unit angle of attack is given by equa- 
tions (6.3) and (6.6): 


C T (a) C-p 4* C T (o ?) . 

Finally the moment coefficient becomes 


Cjy[(a) = (Cyj;^ + C Mt ) sin a cos 


(6.7) 


where Cy[ T is given by equation (6.5). 
is given by 


The neutral point of the force 
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The Numerical Calculation Procedure 


The excessive amount of calculations required for the solution of 
the problem necessitates the use of fast electronic computers. The 
problem was therefore coded for the CDC 3200 and the UNIVAC 1108. 

Before the problem is discussed in detail, the preparation of the 
necessary input of body coordinates and elliptic integral tables is 
presented . 

For the integration of the integral equations, we employ the 
quadrature procedure of Gauss [ 19 ]. The abscissas of the quadrature 
are the zeros of Legendre’s polynomials. The body under investigation 
as sketched in Figure 1 is subdivided into NT characteristic sections. 
The JT-points of the quadrature procedures are distributed within each 
section; in Figure 1, we set JT = 8. In this manner NT X JT points are 
placed along Fx 5 Fa and By choosing the right distribution of 

characteristic sections, it is possible to place a larger number of 
points in the vicinity of the separation point. For instance, in 
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Figure 1, sections 2 and 3 are smaller than the others. The separation 
point was always an end point and starting point of a new section. For 
an abrupt separation, this section remained constant; however, for 
smooth separation, the point s^ is for a constant base pressure coef- 
ficient obtained by an iteration. 

The complete elliptic integrals G n (k 2 ) and F n (k 2 ) are calculated 
with the aid of Landen*s transformation Tit]. This avenue of calcula- 
tion proved to be the shortest and best. Only a few iterations were 
necessary to obtain an accuracy up to 12 digits. In most cases the 
number of iterations remained below 8. 

The geometry of the forebody and the assumed cavity are read into 
the computer in terms of the characteristic sections s^ and r^, where 
s^ and r^ are the arc length and the radius of the particular charac- 
teristic section, respectively. A subroutine specifically developed 
for the particular body distributes the points s v of Gauss* quadrature 
procedure within the characteristic sections up to the separation point. 
In Figure 1, there are, for instance, two characteristic sections, 1 
and 2, with three boundary values and i = 16 points; 8 for each 
section. For these points the subroutine determines the. values x v , r v , 
Xyj x", r!J,, K y and the angle of attack a v = 1 . For the cavity or 

wake, starting with the point j = 17 in Figure 1 and ending with N = 48, 
a table with (N-i) values is read into the computer; these are s^, r^*, 
OCy % With these initial data, the following are calculated: r*, xj,*, 

r"*, x*, x"*, K* . These data are either assumed or approximated. 

We can now start to determine the source distribution q Q (s) of 
the axial flow case. 

(a) We consider r% as a streamline on which condition (1.21) is 
satisfied. With equation (3.1) and n = 0, we obtain the set of N 
linear algebraic equations which is solved for (q Q )* and which satis- 
fies the boundary condition (1.21) in N points on Ti + Fg. 

(b) The next step leads to the determination of the velocities 
( u o)v> (v*) v with the help of equations (3.3) and (3.4), where we set 
n = 0. Thereafter, we calculate (Cp) v and the perturbation potential 
(cp*) v* We add to the perturbation potential the corrective terms of 
equation (4.6). 

(c) Then we solve the system (4.11) with the corresponding right- 
hand sides twice. We obtain the complete solution (4.12) with an 
unspecified constant C. The partial solution (Z^q Q ) v is normalized to 

C = 1. 
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(d) The function (q Q ) v is inserted into equation (4.16), and the 
constant C is changed in such a fashion that the tangential velocity 
Vxj at v = j downstream of the separation point sg is equal to the 
specified tangential velocity V^B • l n Figure 1, j was set equal to 
17. If the condition for V-pg is satisfied in one point on 12 then it 
is satisfied approximately on the entire length of r 2 * 

(e) The general solution q ov now with specified constant C is 
inserted into equation (4.15), and a normal velocity Vn v is calculated 
on r*. With the equations (4.19) and (4.20), a new line r** will be 
calculated. This line is then the new free streamline for a new itera- 
tion step, and we start at point (a). If the absolute value of the 
difference of |Cp(s) ~ Cpg | is equal to or less than .001, we terminate 
the iteration process and consider this case as solved. 

For the body under normal flow, similar conditions prevail. A 
rough approximation of the angle of attack distribution was already 
stored in the machine. The geometry of the cavity stays fixed as 
obtained in the zero angle of attack case; therefore, the Kernel of 
the integral equation does not change during the angle of attack 
iteration. 

(a) We consider for the first step of the angle-of -attack itera- 
tion the cavity or wake as a body with solid contours. This makes f* 
a streamline. The normal velocities diminish on the contour r*- We 
solve for the source distribution p.* which satisfies, for the given 
angle-of -attack distribution, the condition (1.21) in N points on the 
body and cavity. 

(b) With this source distribution jj*, we can calculate the veloc- 
ities u* v , v* v , and w*^ v with the aid of equations (3.3) through (3.5), 
the local normal force C^y with equation (5.22), and the perturbation 
potential with equation (5.29), where cp*J is given by equation (2.2). 

(c) The solution of (5.32) provides us with the complete expres- 
sion of the source strength with a partial solution Ap. v based on C = 1 . 
The partial solution A| i Q needs only to be calculated once, since neither 
the right side nor the elements of the matrix change any more throughout 
the iteration. 

(d) We determine the constant C which makes the tangential veloc- 
ity Vxxb disappear on F* at the point y = j , equation (5.35). Thus, 
the tangential velocity will approximately disappear over the entire 
length of F 2 j since" the potential function has changed accordingly. 
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(e) The constant C is determined, and one inserts the source dis- 
tribution jj, into equation (5.36) and calculates a normal velocity com- 
ponent V^ 1 along the cavity. Finally, one can compute an improved 
angle-of-attack distribution along r% which represents the basis for 
a new iteration step, e.g., (5.39). If the absolute value of the 

normal force coefficient is |c*(s,a)| ^ .001, we terminate the pro- 
cedure and the problem is solved. 

This represents the basic procedure of the iteration scheme. In 
those cases where other models are employed, the corresponding equa- 
tions have to be used; otherwise, the basic structure of the procedure 
remains the same. 


IV. CALCULATED CASES AND COMPARISON WITH EXPERIMENTS 


In order to gain an insight into the magnitude of the error which 
is involved by neglecting the infinitely long wake or cavity and by 
the numerical procedure, the flow field about a half-body with the 
potential j2f(x,p*,cp) = x-E/ (4jrp*) is calculated, where E is the strength 
of the source at x Q on the x-axis. The strength of the source is deter- 
mined in such a fashion that the radius of the body at infinity is equal 
to unity; then the distance of x Q from the stagnation point on the x 
axis is equal to one-half the radius. The spherical coordinates, the 
radius p* is given by p* = sin (cp/2)/sin cp. The pressure distribution 
is obtained by a simple calculation: 


P-Poo 

C = = 1-4 sin 2 (cp/2) + 3 sin 4 (cp/2). (V.l) 

P 9oo 

The contour of the body is thus given by 


= p* sin cp; 


x o " 


cos cp. 


(V.2) 


Since the numerical method employs the arc length s as the independent 
variable, we would like to represent the half -body in terms of the arc 
length: 


s 



n/~P* 2 + (dp'Vdcp) 2 dcp. 
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We obtain finally, as the arc length, 


s 


sin (cp/ 2) r 
cos (cp) V 


3 . 

4 S] - n 


2 2 


2 


+ F(f;k) - E(|,k) , 


(V.3) 


where k 2 = 3/4 is the modulus of the incomplete elliptic integrals 
F(k,0) and E(k,0). The integrals are known and tabulated in reference 
19. For comparison with the numerical method, we need the first 
derivatives 


x 


f 


£P = 

xp 


cos cp + Si ~2 ^ tg (cp/ 2) 
J l+T- tg 2 (cp/2> 



sin cp - C °2 ^ tg (cp/ 2) 
/l + ^ tg 2 (cp/2) 


(V.4) 


and the second derivative: 


d 2 r 


ds2 p*«VT + (tg 2 (cp/ 2)/4) l 4 cos2 ( c p/2) 2 


sin cp , cos cp , / 

+ ““ 9 ^(9/2) - sm 
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tg (cp/ 2 ) (cos cp + ^ y tg (cp/ 2 ) 


SB 


\ 


8 cos 2 (cp/ 2) ( 1 + 4 t g 2 (cp/2) 


(V.3) 


After about 20R the body was truncated. The bow of the body was con- 
sidered to be solid, r 1? and its separation point was placed at Sg = 1. 
The contour of the half-body downstream of the separation point was 
considered to be the free streamline The pressure distribution was 

given along r 2 an d was identical with that of the half-body as given by 
equation (V.l). If we insert now the exact geometric functions, equa- 
tions (V.2) through (V.5) of the free streamline f 2 into the numerical 
method, and if we compare the two pressure distributions as obtained 
by the numerical method and equation (V.l), we should obtain the error 
of the method. The stations of the characteristic sections s a( s ), 
together with the pressure coefficient at the separation point Cpb(s-d) = 
-.29422, are listed in table 1. Following the arc length s of the 64 
points along r x and r 2 are t1ae source strength q Q (s), the error (ERR) 
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(CpB ( s ) “ Cp(s)/Cpp(s) on t ^ Le gi ven pressure coefficient CPW = Cp B (s) 

of (V.l), the remaining normal velocity VN = V^(s), the calculated pres- 
sure coefficient CP = c|5(s), the perturbation potential cp Q (s) = PO, and 
finally the integral of equation (4.8). The error ERR becomes larger 
and larger toward the end of the truncated half-body, as we expected. 

The error, however, is even in the last point on r 2 > less than 1.63 per- 
cent. The normal velocity VN is almost negligible. In the second part 
of table 1, the geometric functions of the exact half body, equations 
(V.2) through (V.5) are compared with those obtained by the numerical 
method. We notice that the radius of the body as calculated with the 
numerical method becomes greater than 1 close to the end of the body. 

The difference, however, is less than AR =.0006R, at the end of the 
body. Looking at the source strength Q = q (s), we notice that the 
strength attenuates rather fast on the nearly cylindrical part of the 
half body. We approximate the decaying curve by C/x 11 , where n > 8. 

The absolute value of the source strength is very small |q (s)| < 3.10" 5 
in the vicinity of the end of the body. Consequently, we surmise that 
the influence of the truncated part of the wake is negligibly small. 

An error estimation of appendix D also shows this clearly. 


The exact drag coefficient of the forebody, which is in polar 
coordinates and according to (6.1), 


9b 

c d°m = 2lt f c p dr (q>) 


r = JL 
B C PB 2 


s in 2 (cp/ 2) - 2 sin 4 (cp/2) 


+ 6 sin 6 (cp/2) 


- r 2 n C 


PB’ 


was calculated for the wake pressure Cpg = -.2942, a radius rg = .7433 
and an angle of cpg = 1 . 67 6 in radians as Cp = .4944. On the other hand, 
the drag coefficient obtained with the numerical method was Cp = .4946. 

The separation of Sg was smooth, and the curvature of the streamline F 2 
was equal to the curvature of the forebody. 

The sequence of convergence for successive iteration steps is shown 
in Figure 2. Plotted are the coordinates of V % for the axisymmetric flow 
and the position of the wake centerlines z x /R for the angle of attack 
case a = 10 degrees. Even with a rough zeroth approximation (a cylinder), 
only about six iterations were necessary to determine the free streamline 
up to about 15 radii behind the separation point Sg. For the angle of 
attack case, the angle of attack distribution changes barely at all 
after nine iterations. The zeroth approximation was tentatively assumed 
to be or (s )/a = 0. 
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In reference 13, Garabedian determined the flow field around a disk 
using the Riabouchinsky model for two base pressures. For the disk with 
C PB = 0 (the infinite streamline model), he obtained the drag coefficient 
Cd( 0) = .827. The corresponding drag coefficient calculated with the 
present method for the same base pressure yielded .824 < C D (0) < .829 
depending on the position of Sg. (The number in parentheses designates 
from now on the base pressure at which the drag coefficient was^ deter- 
mined. ) For the calculation, the boundary condition was satisfied at 
discrete points along the arc length of the contour; 16 points were dis- 
tributed on the forebody and 64 points on the free streamline up to a 
downstream distance of 60 radii. The points were more closely placed 
at the separation point; 8 points were placed upstream and 8 points 
downstream of Sg. The maximum error of the remaining normal velocity 
was on the order of Vj^j < .001 up to about 30R behind the disk. Further 
downstream, the error in the normal velocity increased to about v N < 

.009. The corresponding discrepancy in the pressure coefficient dis- 
tribution amounted to about |Cp(s) - Cpg | < .001 up to s « 30R. In 
table 2 the calculated pressure distribution and the geometric func- 
tions for the cavity behind a disk for Cpg = 0 are tabulated. We 
notice that in the region of 30-60 radii the curvature frequently 
changes the sign; the value, however, stays relatively small so that 
the influence of the error on the pressure distribution of the forebody 
remains small. Similar fluctuations occurred in the source strength 
q Q (s) in the same region. 

The pressure distribution of a disk approached normally by the 
flow is plotted in Figure 3. At the stagnation point, the pressure 
distribution starts with a horizontal tangent. At the separation 
point Sg, the pressure coefficient has a vertical tangent on fi; on 
the free streamline, however, it has a horizontal tangent. The 
derivative of the pressure coefficient Cpg(s) along has a square- 
root singularity similar to the curvature on r 2 . 

Figure 4a shows the difference in the location of the free stream- 
lines of plane and rotational symmetric flow for zero base pressure, 
and shows one main feature of the rotational symmetric flow field. The 
wake is narrower than that of the plane flow field. For different 
truncation positions Sg, the free streamline f 2 was calculated. The 
marked points are those which were used in the numerical calculations; 
we hardly notice any influence on the position of f 2 at least up to 3 r 
behind the separation point. 

An asymptotic representation of the curvature K(s) and the second 
derivative of the streamline radius r"(s) of the axisymmetric as well 
as the plane case is shown in Figure 4b. The two asymptotes were 
determined by equation (4.26) and (4.29). The constant a 1 was cal- 
culated to be a-L = 4.526. Since also the following points closely 
follow the asymptotes, we can surmise that the assumption of the square- 
root singularity and the representative expression for the asymptotes are 
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justified. The second derivative of the radius of the free streamline 
r 2 approaches, for vanishing e = (s-s^), the value -a-j/2 = -2.26. The 
circled points are calculated. The first point is placed about 
e = 4.10“’ 4 R downstream of Sg. With the help of this point and the one 
following in the r (s ) -distribution along r 2 , the constants a-L and b x 
are determined. The curvature and the other calculated functions were 
not always so smooth as in Figure 4b; however, if the upstream distribu- 
tion of points was made approximately the same as the downstream dis- 
tribution, the functions were always smooth. The reason for this is to 
be sought in the determination of H 0 ( s g) hy extrapolation of the upstream 
values of the strength. 

The second drag coefficient, obtained by Garabedian with a different 
method, was for a base pressure of Cpg = -.2235. The modified drag coef- 
ficient Cp(-.2235) = Cjy/ (1-Cpg) had the bounds .85 < Cg < . 88 . His best 
estimate was Cg = .865. A recalculation of this case produced Cg = .830, 
a rather marked difference, though the maximum cavity diameter in both 
cases agreed much better: D^/ 2R = 2.3 as compared with 2.28 obtained by 

the present method. Figure 5a shows the shape of the forward half of 
the cavity, and Figure 5b shows the pressure distribution over the front 
part of the disk and a part of the free streamline r 2 along which the 
pressure was assumed constant. For our calculation 16 points were dis- 
tributed on the meridian of the body Pi and 40 more on the forward half 
of the free streamline. An approximation of the position of the sym- 
metry point s R was obtained with the aid of the dissipation model. The 
point Sg, where r 2 joins r 3 , was used as the zeroth approximation. With 
the procedure discribed in chapter 4 (section (4)) the calculation was 
continued. The half-length of the cavity was L/Dq 6.4. 

Figure 6 compares the drag coefficient obtained from water tunnels 
[ 16 ] with that obtained by the present numerical method. Cones of dif- 
ferent vertex angles were placed into a water tunnel, and their bases 
were ventilated to obtain different base pressures. The theoretical 
results reproduce perfectly the experimental data. Some numerical 
results for the drag of cones are given in table 3. The numbers 
designated with an asterisk are obtained with the Riabouchinsky model; 
whereas, the remaining ones were obtained by the dissipation model. 

All models, despite their differences, converge to the infinite stream- 
line model for a base pressure of Cpg =0. A comparison of the two drag 
coefficients for Cpg - -.41 and the disk, (3 = 90 degrees, gives an indi- 
cation of the difference between the Riabouchinsky and the dissipation 
model. However, for smaller Cpg, the difference in the drag coefficient 
will definitely become greater. 

Figure 7 shows the pressure variation in the wake of a disk as it 
was measured in the water tunnel [ 20] and in the wind tunnel [ 21 ]. In 
reference 20 different torpedo head forms were tested. In the case 
shown, a blunt cylinder in axial flow was tested; therefore, the com- 
parison is not quite justified, since the thick cylinder had a tremendous 
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effect on the wake. A better comparison was obtained from reference 22. 
On the whole, both wake pressure distributions have the same character. 

The theoretical curves were obtained with the dissipation model for a 
base pressure of Cpg = -.41. Across the front part of the disk, again 
the pressure agrees very well with the tests. Along the wake, the 
theory is right. only in tendency. The shape of the pressure distribu- 
tion suggests that a stream tube with a smaller radius than that of 
the dissipation model should be devised. For this model, the free 
streamline will be followed a little bit longer until it meets the 
line At this point a free stagnation point is formed. 

Two additional points obtained from wind tunnel tests of 1 reference 
21 were drawn into the drag plot of Figure 7. The Reynolds number of 
these tests was 7 X 10 4 , The drag coefficient obtained with the numeri- 
cal method and the Riabouchinsky model for Cpg = -.41 was Cq = 1.176. 

The drag coefficient obtained with dissipation model and the same Cpg 
was C D = 1.177. 

The cavity behind a disk, photographed in reference 22,, is repro- 
duced in Figure 8. The measured base pressure, Cpg = -.188, was remark- 
ably constant throughout the cavity up to a distance of six disk diameters 
downstream. The free streamline for the same Cpp was drawn on the photo- 
graph, and a very good result was obtained. The free streamline corre- 
sponds to the dissipation model, and therefore the streamlines run 
parallel from the maximum cavity diameter on downstream. 

The flow field and the wake behind a disk with a base pressure of 
CpB = “*41 were surveyed in a wind tunnel [ 21 ]. The pressure and veloc- 
ity profiles in radial direction were measured, and the corresponding 
mean streamline i|r - 0 was calculated. The theoretical free streamlines 
obtained with two different models were compared with the experimental 
one in Figure 9. The measured drag coefficient was given as Cp(-.41) 

= 1.13, whereas the theoretical yielded Cp = 1.18. The free stream- 
line of the Riabouchinsky model was not plotted beyond x/R ^ 3, because 
the measured pressures along the wake centerline differed greatly from 
the assumption of constant pressure for the theoretical model. 

We now come to the discussion of separation points for convex 
bodies, especially for spheres. For this purpose we refer to the 
schematic representation of possible separation points in Figure 10. 

On the surface of the sphere exists one region (S-A) beginning at the 
stagnation point Sin which any of the separating free streamlines of 
the infinite wake (Helmholtz model) are convex and have infinite curva- 
ture at the separation point S. The streamline intersects the body 
contour within S-A, unless the body ends at s^ . The free streamline p 2 
separating smoothly at point A, has a finite curvature, which is equal 
to the curvature of the body at s-g. The maximum wake occurs for a 
smooth separation. The separation angle belonging to the point of 
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smooth separation is designated as 0 S (smooth). Any free streamline 
separating between A and B has a point of inflection. Immediately 
downstream of the separation point r 2 concave. A distance suffic- 

iently downstream of Sg, r 2 will change the sign and become convex along 
the rest of r 2 * The free streamline separating at point B approaches 
the x-axis asymptotically without changing the curvature of r 2 , which 
is concave. The streamline forms a cusp at infinity, and therefore 
has a free streamline width of zero at infinity. 

The free streamlines separating between points B and S' all form 
a cusp on the x-axis a finite distance downstream of the separation 
point. The tangential velocity V^b an <3 T 2 this region decreases 

with increasing separation angle 0 S and becomes zero at the rear stagna- 
tion point S 1 . These cavities have zero drag. The corresponding cavi- 
ties with a finite cavity width have negative drag. 

A selection of different separation angles and the peculiarities of 
the associated free streamline are given in the following figures, which 
we will discuss mutually. Figure 11 shows the behavior of the free 
streamline r 2 immediately downstream of the separation point sg (smooth), 
0 Q (smooth) = 57.57 degrees. The cycled points are calculated. The 
solid lines represent the asymptote for the curvature k( s) and the 
second derivatives x rr (s) and r”(s). The curvature approaches k q = 

K(sg) ~ “l* The other two functions are consequently finite and 
approach their respective x^ and rj^ 'value of the contour p 1 . The 
infinite wake from a point s^ = 50R was neglected. The procedure to 
determine the smooth separation point was as follows: The character- 

istic section downstream of the separation point sg was for two separa- 
tion angles 0 Q - 57 degrees and 0 S ~ 58 degrees specified to be the same, 
= 0.01; thus, the first point y = j had the same distance (sj - Sg) 
for both cases. The constants a 1 which were obtained were plotted versus 
0 S and interpolated linearly for that 0 S for which a 1 = 0. This interpo- 
lated value was used for the next approximation. For a sufficiently 
small value | a - l | < 10“ 9 , a x was arbitrarily set equal to zero and the 
corresponding angle 0 Q was defined as 0 S (smooth). The value of the 
constant b x was then b 1 ~ .8817. In addition to the geometric func- 
tions the source strength q Q (s) is plotted. The constant c-j_ vanishes 
correspondingly. The asymptotic expressions for the geometric func- 
tions, as well as for the source strength, represent the corresponding 
numerically obtained functions very well. 

Figure 12 shows the indeterminancy of the separation point on a 
sphere for a base pressure coefficient of Cpb = 0* Three separation 
points are chosen. The corresponding pressure distributions on sphere 
surface and free streamline in the neighborhood of the separation points 
are plotted in figure 12a. The radial derivative with respect to the 
arc length is plotted in figure 12b for the three cases, and figure 12c 
shows the corresponding curvature k(s) of the free streamlines. The 


83 



curve with the index 2 indicates the case for smooth separation. Here 
the derivative of the pressure coefficient with respect to the arc 
length dCp/ds is a continuous function on the sphere and free stream- 
line, the condition for smooth separation. At this point the radial 
derivative and the curvature are smooth functions of the arc length. 

The curve with index 1 shows the case where the flow separates upstream 
of the point of smooth separation. The pressure coefficient Cp(s) on 
the sphere vanishes with an infinite negative tangent at the separation 
point, and the curvature k is convex and greater than the curvature of 
the sphere. Consequently, the free streamline penetrates the solid 
body, and this case of an abrupt separation is therefore physically 
unrealistic. One other abrupt separation is shown by the curve with 
the index 3 downstream of the point of smooth separation. The pres- 
sure vanishes at the separation point with a positive infinite tangent. 

The minimum pressure occurs on the sphere, and the pressure in the cavity 
is higher. The curvature of the free streamline at the separation point 
is concave for a short distance, has an inflection point, becomes convex 
further downstream, and decreases monotonically . For cavitational flow 
one can argue that the pressure must be a minimum in the cavity because, 
otherwise, even a small reduction in the pressure coefficient Cp would 
induce cavitation elsewhere. This implies that the free streamline must 
be convex toward the cavity. By Bernoulli 1 s theorem, this is equivalent 
to the condition that the velocity is a maximum on the free streamline. 
Assuming convexity, for the free streamline not to penetrate the obstacle, 
the free streamline must have finite curvature. In fact, the local curva- 
ture of the obstacle cannot be exceeded. These Brillouin separation 
conditions point out that smooth separation should be the only physical 
possibility for cavitation to occur. 

The curvature k( s) of some selected free streamline P 2 is shown in 
figure 13. For separation angles from 0 S (smooth) to « 110 degrees, 
the base pressure coefficient was arbitrarily selected as Cpg = 0. For 
higher separation angles 0 S > 0^ the pressure coefficient was limited 
and could no longer be chosen freely. We notice that the free stream- 
line has a point of inflection within the region A-B ( 0 S (smooth) < 0s 
< 0 S ) which shifts downstream on with increasing 0 S . The concave 
region increases thus until at 0^ « 110 and Cpg =0.0 the free stream- 
line is only concave and is therefore forming a cusp at infinity. If 
the separation point shifts beyond B, the corresponding free streamline 
is only concave. In Figure 14a, we notice that the wake or cavity width 
decreases with increasing 0 S . For 0 S > 0^ the free streamline r 2 anc * 
the streamline P 3 of the dissipation model have a point of contact of 
first order. The diameter of the stream tube T 3 vanishes for a certain 
value of CpB- We designate this value to be Cpg (lim). The free 
streamlines then form a cusp on the x-axis, and wake form I is obtained. 
The particular (Cp^, 0 S ) -combination is plotted in Figure 14b. For 
CpB values below the line designated as M Cpg (lim)" the free stream- 
lines intersect the x-axis and each other. On the other hand, if we 
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choose C PB values above Cpg(lim), we obtain the wake form II with nega- 
tive drag. On the line Cp B (lim) the drag coefficient is zero because 
the corresponding wake width is zero. For 0 Q > the free streamline 

P 2 is only concave, and one can prove with equation (C.16) that Vp B < 1. 
Therefore, only positive pressure coefficients are possible in this 
region. Since the mean curvature increases with increasing 0 S , the 
tangential velocity Vp B on P 2 decreases until it becomes zero for 
0s = 180 degrees. Figure 15a shows, in addition to wake form I, also 
wake form II for the case 0 S = 130 degrees. The pressure coefficient 
Cpp(lim) = .222 produces approximately zero drag. However, if we 
choose Cp B = .25, we obtain a finite wake width with the corresponding 
drag coefficient of = -.026. In order to obtain the cusped wake 
form, the rearward region of the cavity had to be changed. The pres- 
sure coefficient on the fixed cone surface (its geneatrix is an expo- 
nential function) did change from the specified Cpp only slightly. 

Only in the region of the rearward stagnation point did the pressure 
change appreciably; therefore, the point of intersection between r? and 
the x-axis is quantitatively and qualitatively not quite right. One 
other possibility, however, to obtain the cusped cavity is the extrapo- 
lation of the wake form II to lower Cp B values. 

Figure 15b shows an additional region where the separation angle 
0 S and the wake pressure cannot be chosen arbitrarily. The correspond- 
ing wake forms for a Cp B = -.1 are plotted in Figure 16. If we hold, 
for instance, the base pressure constant (Cp B = -.1) and let the 
separation angle 0 S change to higher angles, then the wake width 
decreases and with it the drag coefficient. The curvature of P 2 is 
concave and becomes convex shortly before it joins the streamline r 2 . 
The convex region finally disappears, and P 2 is only concave. A con- 
cave streamline cannot sustain Vp B > 1, and therefore the streamline 
model breaks down. For Cpg = -.1, the separation angle 0 S ^99.6 
represents the limiting rearward angle possible for the model. One 
other try for 0 S = 99.75 failed. Figure 13c also shows that, for 
decreasing 0 S , the contact point s^ between an d T 3 moves forward. 
With decreasing Cp B , the separation angle 0 S (smooth) moves to higher 
angles. The limiting separation angle j^*, however, moves forward from 
110 degrees for Cp B = 0 to a position of 90 degrees for a Cp B ~ -.525. 
At this point, the free streamline has disappeared and p ± leads 
directly into p 3 . These investigations could be carried on almost 
without limit. Although there are still a number of interesting ques- 
tions about separation from round bodies, we are presently content 
only with this limited investigation. 

In Figure 17 an attempt was made to predict the laminar separation 
point of the sphere. (The boundary layer method used can be found in 
reference 23.) The boundary layer equation is solved by a series 
approximation of the "separated" velocity distribution, by the use of 
an iterative procedure. For a fixed base pressure coefficient of 
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Cpb = -.45, an arbitrary separation point around 0 & 80 degrees was 

chosen. The resulting velocity distribution on the sphere was developed 
into a polynomial containing terms of up to s 7 , where s is again the arc 
length. The boundary layer equations were then solved, and a new separa- 
tion point was obtained. The wind tunnel experiment yielded a separation 
point of 0 S —81 degrees; the iteration procedure predicted a separation 
point of 0 S — 77 degrees plus. One of the difficulties was to reproduce 
the velocity distribution in the neighborhood 6f the separation point. 

A more modern procedure might have given a much closer result. Fig- 
ure 17 shows the theoretical pressure distribution for a laminar separa- 
tion at 0 S = 80 degrees, and for turbulent separation at 0 S = 143 degrees. 
Due to the laminar separation bubble (T) on the lee side of the sphere, 
the theory does not quite reproduce the actual pressure distribution 
over that portion of the surface. The theoretical pressure distribution 
was selected in such a fashion that the pressure distribution of test 
and theory coincided very well over the forward part of the sphere. 

For a separation angle of 0 S ^ 130 degrees, we obtain the lowest pos- 
sible pressure coefficient as Cpg(lim) = .222 with Cp = 0. 

The Brillouin separation conditions were tested in a water tunnel 
[24] on a sphere, (The results are plotted in Figure 18.) The drag 
coefficients for spheres obtained by the use of axial singularities 
(flow past a half-body) were given as Cg(0) = .30 in reference 24. This 
result can be compared to the smooth separation occurring at an angle 
of about 0 S = 57 to 58 degrees measured from the forward stagnation 
point. The corresponding drag coefficient obtained by the present 
method yields C D (0) = .31. It seems, however, that cavitation does not 
occur at the point of smooth separation (Brillouin point), probably 
because of surface tension. The cavitation occurs, rather, at greater 
angles 0 S ^ 80 degrees depending mildly on base pressure coefficient. 
Using these separation angles, the corresponding drag coefficients were 
calculated. Varying the separation angle by a few degrees and plotting 
the corresponding drag result, it was found that the drag coefficients, 
shown in figure 18a, were closer to the minimum drag coefficient of the 
sphere for fixed base pressure coefficient. 


Figure 19 shows the behavior of the free streamline immediately 
downstream of the separation point Sg = 0 Q = jt/ 2 from a sphere. For 
this point r 1 = 0, r M = -1 and xj = 1. The curvature of the sphere is 
constant and K Q = - 1. With equation (4.29) and only a x considered, 
we obtain, for e 0, 


K(S) 
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where a ± = .3616. This function is plotted in figure 19 as an asymptote. 
Close to the separation point the free streamline is very well repre- 
sented by the asymptotic expression. For this point, only the second 
derivative of the radius r(s) contributes to the square root singular- 
ity. Consequently, the function x M (s) must approach a finite value: 

lim x" (s ) = - —k' ~ -.1808, 

e-» 0 

which indeed will be approached. 

The lift of cones was measured in the water tunnel [25]. Fig- 
ure 20 plots the lift coefficient gradient dCq/da for zero angle of 
attack and a base pressure of Cpg = .0 and -.1. The base pressure, a 
function of the angle of attack, was held constant for these investiga- 
tions. One notices that, with decreasing base pressure coefficient, 
the lift coefficient gradient also decreases. For very small cone 
half-angles, the lift slope approaches the "slender body" result of 
dC-^/da = 2. With increasing cone half angle, however, the lift grad- 
ient decreases and becomes negative until at p = 90 degrees only the 
drag component contributes to the lift. 

In addition to the theoretical lift coefficient gradient obtained 
by the numerical procedure, some approximate theories are proffered 
which were obtained from reference 26. The curve for the modified 
slender body was obtained in the following way. For calculating the 
normal force of a slender body, only the change of the momentum in 
planes normal to the body is considered. For slender bodies, this 
plane moves with the velocity of along the body axis , and the 

result is the well known Cjjx = 2.0. For cones with an appreciable 
half-angle |3, the velocity along the cone surface will be less than 
Uoq. We assume that this velocity will be Vp cos p. We obtain with 
section 6, for the lift gradient, 


C L a “ C % - C D 0 » 


where 


C Nq; = 2V^ cos 2 p = 2(1 - Cp o ) cos 2 p. 


Finally, we obtain with 


% = 2(1 - Cp o ) cos 2 p - C Dq 
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an expression which shows a relatively good coincidence with tests for 
all half-angles including the disk. 

One other method, termed "strip theory," was taken from reference 
26 and plotted into Figure 20. The theory resembles a correspondence 
theory, where two-dimensional results are applied to bodies of revolu- 
tion. 


The lift, moment, and drag coefficient of cones with half-angles 
of p = 15, 45, and 90 degrees and a constant base pressure coefficient 
of Cpg = -.1 are plotted as a function of the angle of attack in Fig- 
ure 21. The disk is here considered as a cone of 90 degrees half- 
angle. For technical reasons the angle of attack was chosen to be 
negative. The drag coefficient C D (a) as a function of the angle of 
attack is fairly well reproduced by the linearized theory except for 
the cone of (3 = 15 degrees, where the theory does not quite predict 
the increase in drag with angle of attack. It is possible that some 
higher terms in a which are neglected here play a certain role. The 
lift coefficient is satisfactorily reproduced as is the moment coef- 
ficient over the angle-of -attack range where separation of the flow 
from the lee side of the body is not expected, a, ^ p. The reference 
point for the moment coefficient was in all cases the center of the 
cone bases and the stagnation point of the disk. 

Figure 22 shows the lift and drag coefficient of a 15-degree cone 
plotted versus the base pressure coefficient. Even for an angle of 
attack of OC = ±20°, the theory is still predicting the right lift coef- 
ficient. Also, the tendency with varying base pressure coefficient is 
well reproduced. At an angle of attack of CC = 0 degrees, the cavity 
was subjected to a hydrostatic lift which was noticed again on the fore- 
body. The theoretical result was therefore shifted toward the experi- 
mental for a, = 0. The change of the drag coefficient with the angle of 

attack is very good for small angles; for higher values of ot, however, 
the theory yields smaller values. This was already noticed in Figure 21. 

In Figures 23 and 24 the same coefficients are plotted as in Fig- 
ure 22 for a 45-degree cone and a disk, respectively. In these cases, 

the cavity had little or no effect on the a = 0 shift for the lift 

coefficient, since the normal force for a 45-degree cone is very small 
and that of a disk is zero. 

The downwash angles behind cones of p = 15, 45, and 90 degree half- 
angles are compared in figure 25. With decreasing cone half-angle, 
which corresponds to increasing lift coefficient of the forebody, the 
ratio of the local wake angle of attack a(s) to the angle of attack of 
the forebody a increases. For a disk, therefore, the wake leaves the 
base essentially in the same direction as the free stream; whereas, 
for cones the downwash angle increases behind the base and approaches 
asymptotically a certain limiting value at infinity, according to the 
normal force of the forebody. 
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To obtain an idea of how good the approximate theory describes 
the wake of a body at zero and small angles of attack, the free stream- 
lines were drawn into the shadowgraphs of the flow configurations of 
Figures 15 and 16. Figure 26 shows a cone of 50 degrees half-angle at 
an angle of attack of a = 10 and 0 degrees. The base pressure coef- 
ficient was measured at Cpg = -.4. The lower edge of the wake is 
fairly well reproduced close to the body base where the shear layer 
of the wake is relatively thin. Figure 27 shows the wake of a 15-degree 
cone at a = 10 and 0 degrees with a base pressure of approximately 
CpB “ -.32. Wake boundary and free streamline show again a fair coinci- 
dence, at least up to two diameters behind the base of the cone. The 
free streamlines for both cones were obtained with the dissipation 
model . 


V. FURTHER APPLICATIONS 


In the area of steady discontinuous flows, other applications of 
the theory are possible. In some engineering problems, the added or 
virtual masses of bodies at separated flow conditions -- for instance, 
meteorological balloons and parachute canopies -- must be known to 
determine their dynamical behavior. The theory yields these quanti- 
ties readily. Furthermore, rotational symmetric bodies can be con- 
structed whose surfaces have a certain prescribed pressure distribution. 
However, the given distribution must obey certain rulesiin order not to 
produce negative body radii. Into the same category falls the problem 
of designing the optimum shape of the cowl of a rotational symmetric jet 
intake. Another class of problems is the flow out of an orifice, jet 
penetration, and cavity and jet flows under the influence of gravity 
fields. The theory is directly applicable to some of the mentioned 
problems; whereas, for certain others, slight modifications in the 
boundary conditions have to be made. 


VI. CONCLUSIONS 


The flow about bodies of revolution with blunt bases separates and 
forms a free streamline which divides the flow field into an outer flow 
region and the wake region with zero velocities. In case of cavita- 
tional flow, this free streamline is clearly visible in flow pictures. 
For wake flows, the free streamline concept is merely an idealization 
of the free shear layer. 
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An integral equation method using singularities on the surface of 
the body and the free streamline was applied to obtain solutions for 
the mixed boundary-value problem. The problem contains two parameters 
for the case of separation from a smooth body: the base pressure Pg 

and the separation point s-g on the surface of the rigid body. With 
the help of boundary layer theory the second parameter can be eliminated 
as was proven with an example (laminar separation from a sphere). The 
base pressure for wake flow must be obtained from experiment. For 
cavitational flow, P-g is equal to the vapor pressure of the liquid. 

An approximate procedure for the body at small angles of attack 
analogous to the zero angle of attack case was developed in the second 
part of the report. The geometric shape of the cavity as obtained in 
the axial flow case was used for the angle of attack case. The sec- 
tions of the wake are shifted to a position where the local normal 
force on the wake is zero. The potential of the bent wake could be 
reduced to that of a straight one after a series development and 
linearization where the angle of attack distribution varies along 
the boundary of the wake. 

A number of examples are given, and theoretical and experimental 
results were compared. Excellent to good agreement of drag, lift, 
moment, and induced drag coefficients was achieved even at moderate 
angles of attack. The excellent agreement between theory and experi- 
ment proves further that the free streamline concept can be applied to 
wake flows, as long as one is concerned with the forces on the rigid 
body only. 
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TABLE 1. COMPARISON OF EXACT AND APPROXIMATE 
PRESSURE COEFFICIENTS ALONG THE HALF-BODY 
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TABLE 1 (continued) 


The Exact and Approximate Geometric Functions of the Half-Body 



Exact: 

Parallel Flow 4- Source 

Approximate Numerical Procedure: 
(The region s H is defined as 
wake . ) 

s 

X 
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RP 
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RP 
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9.9*786- PI 

1 . 1 3 U 7 6 - U 2 

5 . 40266 - 04 

7. 63316 

00 

7.18766 00 
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7.92146 00 

9. 97O7E-01 

6.22706-04 

2 . 41 3 /E-U4 

9. 05106 

DO | 

8.6056E 00 

9. 90116-01, 
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9.59336 

on i 

9.147B6 00 
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TABLE 2. THE PRESSURE COEFFICIENT AND THE GEOMETRY 
OF THE CAVITY BEHIND A DISK FOR C PB = 0 
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LEGEND: 


- arc length 
= x- coordinate 
= wake radius (r(s)) 
= dr/ds 
= d 2 r/ds 2 


CUR = curvature 

Q = source strength 

VN = remaining normal velocity 
on the streamline 

CP = pressure coefficient 
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TABLE 3. CONE AND SPHERE DRAG COEFFICIENTS 


Cone Drag Coefficients 


C PB 

.0 

-.05 

-.1 

0 = 15° 

.1460 


.2014 

45° 

.5042 


.5694 

63.4° 

. 6679 



VO 

o 

0 

.8238 

. 8664 

.9092 


(Riabouch. 


-.2 

-.2235 

-.3 

-.41 

. 6360 




.9955 

1.0144* 

1.0822 

1.1777 

1.1767* 


Sphere Drag Coefficients 




Mod . ) 

-.6 


.3455 
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TABLE 4. THE APPROXIMATE CALCULATION 
OF THE LAMINAR SEPARATION POINT FROM A SPHERE 


s 

v T 

Calculated 

Error 

.027423 

.055215 

.040701 

.2628 

.140419 

.203833 

.207253 

-.0167 

.327660 

.466285 

.471754 

-.0117 

.563907 

.769968 

.769153 

.0010 

.817262 

1.040140 

1.032990 

.0068 

1.053509 

1.218332 

1.221146 

-.0023 

1.240750 

1.289778 

1.300160 

-.0080 

1.353746 

1.285990 

1.285045 

.0007 

1.382162 

1.274379 

1.269342 

.0042 

1.386253 

1.272036 

1.266597 

. 0042 

1.393031 

1.267721 

1.261764 

.0046 

1.401584 

1.261152 

1.255147 

.0047 

1.410755 

1.252393 

1.247384 

.0039 

1.419308 

1.241234 

1.239498 

.0013 

1.426086 

1.227846 

1.232790 

-.0040 

1.430177 

1.212065 

1.228540 

-.0135 


u7 = 1.4845, = -.4385, = .2094, u y = -.0706 

The determination equation for the separation point is: 

1 - .6832s!; + .4666s* - .2198s° = 0. 
£>£>£> 


Solution: 

S B = 1>342, 0 s = 77 ° - 
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(c) The Bent Wake and Additional Notations 


Figure 1. (Continued) 













(The marked points are calculated) 
o r 2 Truncated after 10 radii 

• r 2 " " 30 " 

A r 2 " " 60 " 

(a) The Free Streamline in the Neighborhood of Sg 

Figure 4. The Free Streamline In Plane and Axisymmetric Flow 
About a Disk for a Base Pressure Coefficient of Cpg = 0 
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Figure 8. The Cavity Behind a Disk for Cpg = -.188 [22] 

(Water Tunnel Test) 



Figure 9. Comparison of the Streamline \|r = 0 Obtained from 

Wind Tunnel Tests and Theory for A Disk with C^- = -.41 

FB 
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Figure 10. Schematic Representation of the 
Separation Points from a Sphere 



Figure 11. The Free Streamline in the Immediate Vicinity of 
the Smooth Separation Point of A Sphere for Crm = 


0 


.2 



(a) Pressure Distribution 



(b) Derivative of the Radius with Respect to s 



s/R 


(c) Curvature of the Free Streamline 

Figure 12. Representation of the Possibilities of 
Separation Close to the Point of Smooth Separation 
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1.07 
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105 

104 

103 

R 102 
1.01 

100 

.99 

.99 

.97 


Figure 16. 









s = (s - Sg) 


Figure 19. The Asymptotic Representation of the Curvature and the Second Derivative of the 
X-Coordinate of the Free Streamline in the Immediate Vicinity of the Separation 
Point of A Sphere 0 = %j2 
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Figure 20. The Lift Gradient of Cones as A Function of the 
Half-Angle p and the Base Pressure C™ 
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Figure 21. Drag, Lift and Moment Coefficients for Different Cones 
As A Function of the Angle of Attack for C pB = -.1 [ 25 ] 
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a(s) .5 



(a) The Dimens ionle 


0=l5 o v 


( X 

a =lo 0 A-v ✓ «° 


Figure 25. The Theoretical 
and the Shifted 








(b) The Axisymmetric Flow About the Cone 


Figure 27. The Shadowgraph of the Flow About A 15-Degree Cone 

with A Superimposed Streamline for C = -.32 

PB 
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APPENDIX A 

THE EVALUATION OF THE ELLIPTIC INTEGRALS 


In the evaluation of the surface integrals, we often encounter the 
integral expressions 


F (k 2 ) 
n v y 


*/2 



cos (2ncp) 
[l-k 2 sin cp] 1 / 2 


dcp, 


(A. 1) 


at 2 

G n (k=) - (-1)“ f 

n J [ l-k 2 s in 


dcp. 


(A. 2) 


For the nominator we choose the expansion 


2n 


cos (2ncp) = cos cp - ^ 2 y s ^ n 9 cos 


2n 


. 2 2(n-i) 


2 n 


4 2 (n-2) 


cp + ( 4 Jsin^cp cos 


9 


(A. 3) 

which, after integration, yields for the expression (A.l) the terms 


F o (k 2 ) = K (k 2 ) - G o (k 2 ) - — [G^k 2 ) + G o (k 2 )] 


(A. 4) 


and 


Fi(k 2 ) = - 1) K(k 2 ) - p E (k 2 ) = - GiCk 2 ) + — [G x (k 2 ) + G Q (k 2 ) ] . 

For the elliptic integral (A. 2), we obtain for n = 0, 1, and 2, finally: 


G Q (k 2 ) = E(k 2 )/k’ . 

1+k g - 2 E (k 2 ) - 2K(k 2 ) 


Gl(k 2 ) - ^ 


k' 


(A. 5) 


2\ _ _L 


G 2 (k 2 ) = 


•|l+k ,2 +12+ 1 ^ k g -j - E (k 2 ) - -|^2k' 2 +6 (1+k 1 2 )+2j* K(k 2 ) 


= 1 - k 2 . 


where k' 
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APPENDIX S 


The Asymptotic Development of the Integrals 

In the vicinity of s -> a, the modulus k 2 approaches unity, and 
the elliptic integrals have for k f 2 g .5 the series expansion 


k 1 63 , * 4 925 t e 


K(k^) 4 " 384 k ' " 7680 k 


+ In (4/k' ) 


1 + 


k'' 


, _± k , 4 + 25_ k - 6 + 

+ 64 256 


(B.l) 


E(k 2 ) = 1 - 


k ' 2 39 , ,4 90 , ,6 

— ~192 k “ 640 k • 


+ In (4/k' ) 


k ' 2 3 2 

JsS — + _ k ' 2 

2 + 16 “ 


+ t5_ k -6 + 

128 • * * 


where k' = 1 - k 2 . 

We are able now to determine the first terms of G n (k 2 ) with 
this series expansion; 


G 0 ( k2 ) 


1 1 - T k ' 2 - k ' 4 - ^ k ' 6 . 

k l- 

1 ,,*,15 


. + In (4/k') (~ k ' 2 


+ 16 k ' + 128 k ' + J 


x(k 2 ) = ^ + l k ’ 2 + k ' 4 +M|i k' S + ... - In (4/k 2 ) (| k ' 2 


+ ii k ,^ + 165 k ,s ^ 

+ 16 + 128 / 


(B. 2) 


(equation (B. 2 ) continued on next page) 
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Gad**) * ^ ' 1 + f k '" + TT k ' 4 + ^ k' S + 

k i- 


576 


- ln(4/k*) (yk |2 + ff 1 k' 4 + yyf— k' 6 + 


In general, one can represent the elliptic integrals for k' 2 g .5 by 
the summation formulas 


K( k=) . i.(w> £ y } 4 + m)*) - ^ (a 

v=l V=1 


J- + .2. 

, 1 • 2 3-4 


(v-1) 1 ,,2v 

(2v-l)2 vj f ( 2 + 11)2 tTtT 
M* 


and 


(B. 3) 


E (k 2 ) = In (4/k 1 ) > 

L I 

v=l 


'/ .| 2 V (V-1) 

^D’.V. *, ' 2 


(■ j + + n)) + 1 - r k |2 


V -1 


v=2 

k- 2 ^~ l) i i 

• (v _i)*. v » n C ~2 + n) (j + n) • 

fi=0 


The previous equations are not suited for the region of very small k 2 . 
The summation formulas for determining K(k 2 ) and E (k 2 ) are therefore 
d if ferent: 


2 2 
•2 + 3-4 + 


+ 


(2v-3)(2v-l) (2v-l)2v 


128 


CO 


K(k 2 ) 



(v-1) i 
n (y + n ) 2 
n=o 


| ' 1 + C-Jc 2 + C 2 k 4 + . . . 


and 


(B.4) 


E(k 2 ) 



+ A 

v=l 


k 2v 

v! v'. 


(v_1) 1 1 

H ("2* + M-) (j + M-) 

H=0 . 


r 

I4<i 1 k 2 +d 2 k 4 + ••• 


After the series have been determined, we insert the coefficients of 
K(k 2 ) and E (k 2 ) into the equation (A. 5) and obtain for the region 
k 2 ^ .5 

r 

k'^GxCk 2 ) = | j (2 [d 2 -C 2 ] - [di- 2 C 3 .Dk 2 + (2[d 3 -C 3 ]-[d 2 -2C 2 ])k 4 + ...J 

(B. 5) 


k' 2 G 2 (k 2 ) 


n 

2 


(16[d 4 -C 4 ]-8[2d3-3C 3 ] 


4- d 2 -8C 2 )k 4 4- 


(16 [d 5 -C 5 ] 


-8 [2d 4 -3c 4 ] + d 3 -8C 3 )k G + ... . 


The Asymptotic Development of the Velocities for Small Cone Angles 


For very small cone angles (p s 15°), the series approxima- 
tions derived in chapter IV. 2, are not very well suited for the deter- 
mination of the velocity components. For instance, the logarithmic 
terms of the velocity components become very large for small r(s); 
whereas, f (s , a) disappears rather rapidly with increasing a. A special 
series approximation was therefore advisable for small cone angles. 

We set again e = cr-s . For a cone, the following relations 
hold: r = r's, p = r*a, and x = x T s, | = x' a. 
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The elliptic modulus k T 2 becomes 


now 


k . 2 = izzlll ± (r-P ) 2 = ef 

(x -|) 2 4- (r+p ) 2 e 2 4-4r 1 2 s (s4e) 5 


(B. 6 ) 


which is inserted into (B.l) and with (A. 5), we obtain finally 


V k > 


e 2 4-4r 1 s (s4-e) 
72 


1 1 1 
4 2 I e 


8 r's 


„ p, € 2 +4r's(s+s) 7 3 , 8 r's 

Gl (k 2 ) » -2 + t - - In -r-T- 


G 2 (k 2 ) 



For simplicity, we retain the argument of the logarithmic terms. 

We now introduce these expressions into the integrand of the 
u n -velocity component and obtain the series 


T 2 r 


f ( s > a >u„ = _c ln( s ) . , 

n L X i/ 2 


1 x 1 _j_ 2 r * x ! s 


eX 


1/2 J 


+ q! <■) 

n X 1 / 2 


(B.7) 


where X = e 2 4 - 4r f s (s4-e) . In this expansion, all terms which were pro- 
portional to e were dropped. The integral with the substitution 


Xg = (s E - s) 2 + 4r ' ssj 


finally becomes 


S E- S 


h f f(s ’ a) u da = - L 2r ' ln 


s -s+ 2 r' 2 s+xi ,/ 2 
i E E | 


2 r ’ s 


(equation (B. 8 ) continued on next page) 
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r"(s +s)-xi /2 -I s 

+ i ( , -sh!-?) i 2r ' s in i- 

E 


s -s-b2r f s+Xi;/ 2 


2r ! s 


Next, we develop the v n -component into a series 


£<s -^ = - 2 r, 2 { q ” <s) fe + i + ^(~ 
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The integration of this expression yields 
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For the circumferential velocity component w n , we obtain 
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and integrating this term, we can write 
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(B. 12) 


With the aid of the preceding expansion for the perturbation velocities, 
we obtain the normal component Vjq(s) as 


f(s,o) v =.4q (s)r 
V N n 


2 

X 


2 2 
S ^ 


[ 1 8r f s 


a+2 "Vj 


(B. 13) 


The integral is 
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(B. 14) 
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Similar expressions for the tangential velocity V-p(s) are 


f (s,a) v = -2r' 2 {q (s) ^-2- + J- + ^ 1+2n _ ln &r^\] + q^s)s'. 
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(B. 15) 

After a closer inspection of these terms , we notice that they are not 
suited for determining the velocities around a blunt cone (disk), since 
r 1 1 and x T 0. 
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APPENDIX C 

THE COORD INATE TRANSFORMATION OF THE VELOC IT IES 


For the tangential and normal velocity components on a streamline, 
we find the following relations: 


V T = (l+u 0 )x' - v Q r' 


(C.l) 


and 

V N = - (l+u Q )r ' + v Q x*. (C. 2) 

Certain relationships exist between the derivatives of the streamline 
and its curvature: 


and 


K = x f r !? - r*x ,f = r ,f /x f = -x n /r 1 
x'x" + r 1 r f! = 0. 


(C. 3) 


(C.4) 


If we differentiate the normal velocity with respect to the arc length 
of a line on which does not disappear, we obtain 


dv du dv 

- — = - — — r ' - (l+u )r" + x' - — 

ds ds v o' ds 


+ v x". 
o 


(C. 5) 


With (C.3) and (C.4) the relation 


(l+u o )r" = (1 +u 0 )x'k; v q x" = -v Q r ’k 


(C. 6) 


is obtained. Inserting, this relation into (C.5) and expressing the 
total differentials by the partial one, 


du 

ds 




(C. 7) 
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and 


ds 


to . ^ts. , 

Sx X Sr r 5 


yields finally, with Su Q /Sr = Sv Q /Sx, 


dV 
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'Su Su 

11(0 o 

— TT 1 V* 1 ' 


ds 


+ V * x r V5r - & 




(C.8) 


Differentiating equation (C.l), we obtain, with the same relations, 
the change of the tangential equation along s: 


dV T o 

T _ ,2 O 

ds x Sx 


Su 


+ 2x' r ' 


Sv 


__! 1 + r ' — 2 . 
Sr + r Sr • 


(C. 9) 


Very simple relations are obtained after a coordinate transformation 
for equations (C.8) and (C.9). A rectangular coordinate system t and 
v has its origin at the separation point s B . The T-coordinate is 
orientated tangentially, and the v-coordinate is normal to the free 
streamline. The normal v is positive if it points to the outside flow 
field (or to the left of the streamline). 

The transformation equation is 


(x-Xg) = x't - r'v; (r-r ) = r'r + x'v. 


(C. 10) 


The derivatives of the velocity potential cp(x,r) with respect to 
the coordinates t and v are 


Sep _ ScpSx ScoSr 

St SxSt SrSx 

and 

Sep _ SepSx SepSr 

Sv SxSv SrSv 


(C.ll) 
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If we replace in equation (C.ll) the partial derivatives by the total 
derivatives of the streamline, 


Sx . 

Sr t 

. 

Sr f 

— = X 1 ; 

— = r 1 ; 

— = -r' ; 

— “ x ! , 

St 

St 

5v 

Sv 


(C. 12) 


we will obtain the partial derivative of the potential with respect to 
the body coordinates x and r: 


and 
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Sx 
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St 


. ^ r t 

Sv 


= ik x » + See r > . 

Sv St 


(C. 13) 


The next higher derivatives are 


& = x' 2 ^ - 2r ' x' + r - 2 & 

; St 2 5vcSt Sv 2 



|f| = r . 2 & + 2r ' x' ^fSL. + x’ 2 ^ . 

ur St 2 StSv Sv 2 


(C.14) 


The change of the normal velocity along the arc length is obtained by 
inserting these expressions into equations (C.8) and (C.9): 




(C. 15) 
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If the general line is a streamline, then the conditions V N = 0 and 
dV^j/ds = 0 hold, and equation (C.15) becomes 



(C. 16) 


The change of the tangential velocity along s becomes an identity: 
dV T /ds = d 2 cp/ds 2 . 

We consider now the tangential and normal velocity of the body in 
normal flow. The tangential velocity is 


V T X = u l xf + (Of + v i) rI 5 
and the normal velocity can be written as 


(C. 17 ) 


V N = -Ujr' + (a + v 1 )x' . (c. 18) 

The total derivative of the normal velocity with respect to s is with 
the relations mentioned above 



The third term on the left side vanishes if the angle of attack a is 
constant along the body. 


138 



APPENDIX D 

THE ERROR INTRODUCED BY NEGLECTING THE INFINITE WAKE 


In all cases in which we employ models with an infinite wake and 
truncate the wake a finite distance downstream of the separation point 
at s E we introduce an error. This error will definitely increase with 
decreasing distance of the end point s E from the separation point s E . 
However, since we specify the pressure distribution along the free 
streamline, the error will not influence the velocity distribution to 
any extent on the forebody and therefore the drag of the body; but it 
will influence the position of the free streamline, the source strength 
q n (s), the perturbation potential, and in the normal flow case, the 
angle of attack distribution. 

The estimation of the magnitude of the error is our next goal. 

Some simplifying assumptions are made for this purpose: 

(a) The wake downstream of the truncation point s £ is cylindrical, 


p = r E = const. 


and therefore, 


a = § + const. 

(b) The source strength along the truncated wake attenuates 
according to a simple law 


q n 00 = q n ( x E ) ( X E/ X ) k (d- 1 ) 

where k is not necessarily an integer. 

The assumption (a) is exact for the dissipation model. For the 
Helmholtz model where V T (s) = 1 on r 2 , (a) is only an approximation. With 

respect to assumption (b), we should mention that Levinson determined 
in reference 17 the asymptotic shape of as 


r 


C* 

(In x) 1 / 4 


r 
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1 In (In x) 
8 In x 


+ 0 (1/ln x) , 


(D. 2) 


139 



where C' v represents a constant which can be brought into a relation 
with the drag of the body. For slender bodies, we can assume that 


/ \ dr 

V x) “h 


n* 


2 \Tx (In x) 1 / 4 L 2 (In x) 1 / 2 


1 - 


+ . • . 


CD. 3) 


The source distribution q Q (x) vanishes with about J A/x. On the other 
hand, we notice, after a close inspection of table 1, that the exponent 
for a cylindrical wake is k ^8. 

In all cases, for the potential as well as for the velocity com- 
ponent, an integral of the form 


J^q(i)k(x, |) d| 
X E 


must be calculated. In order to carry out the integration along a 
finite distance, the integral boundaries have to be subjected to a 
transformation. We therefore set 


T] 


1 d£ - ft 2 

Ji p = _dT l ; 6 




and 

a 2 = (r+r E ) 2 + x 2 , 

X = r) 2 a 2 - 2xr) + 1, 
X E = - 2xr] E + 1. 


(D.4) 


(D.5) 


With these substitutions, we obtain the expression 


(|-x) 2 + (r+r E ) 2 = [T) 2 a 2 - 2x^+1] - ^ X, 


(D.6) 
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and the elliptic modulus becomes 


k 2 = 


4rr n 2 


7] a 2 - 2xt] 4* 1 


= 4rr t] 2 /X. 


(D.7) 


We start with the evaluation of the potential, 
an infinite wake the perturbation potential is 


For the model with 


V S > = “ 2^ J q n (CT) 
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The corresponding expressions for the perturbation velocities are 
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V s > " 
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(D . 1 1 ) 


The additional tangential velocity component AV T = Au n x* + Av n r f is 


AV t (s) = 


q (x^ ) r 
^n v E E E 
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T) k (xr)- 1 ) x ' G n (k 2 )+rj k+1 r ' [rG n (k 2 ) - r E (G n+1 (k 2 ) +G n _ ± (k 2 ) ) / 2 ] 

X 3 / 2 


dq. 


(D . 12) 

The additional normal velocity component AV^ = -At^r 1 4* Av n x* * 
which is induced by the neglected part of the wake, is given as 


AV n (s) 


q n< x E )r E X E 



o 


- q k (XT)- 1 ) r ' G n (k 2 )+T) k+1 x ’ [r G n (k 2 ) - r (G n+1 (k 2 )+G n _ x (k 2 ) ) / 2 ] 

X 3/2 


d > 


(D . 13) 

We replace the elliptical integrals in equations (D.8) through 
(D.13) by their respective series from appendix B and consider only 
the first few significant terms: 


142 



V kE > -f C 1 + T + •••) 
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G o (kE > "I ( L + f kE+ •••) 

Gi(k 2 ) = J (j k2 + . 

G 2 (k 2 ) = | (0(k 4 ) + ...). 


(D. 14) 


If we assume that the expression for X is only a little different than 
unity (X « 1), we can develop the express ion . for X into a series. We 
obtain 


« 1 - \ (a 2 T) 2 ) + XT] + ... 
x l/2 2 

— ~~ ^ 1 " 7 C a2r l 2 ) + 3 x T] + ... • (D.15) 

X 3/ 2 2 

We determine the additional normal velocity AV^j(s) on rp in the vicinity 
of the separation point by inserting the expansions (D.14) and (D.15) 
into (D.13). The most important term under the integral is then 


av n ( S ) 


q n (x E )r E X E 


— rj (r 1 + [x 1 r-r * x ] t]) (1 + 3 xr)) dr\. 


(D . 16 ) 


The integration of this expression and inserting the proper integration 
boundaries yields the additional normal velocity 
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(D. 17) 


AV n (s) 




_k+i 


J. x* r + 2r 1 x 

(k+l) x| 


For very large x^, we can also neglect the second term within the brackets. 


AV n (s) 


, v r E r T 1 
^n^ X E^ 2 k+l x r 


(D. 18) 


The error is thus proportional to the source strength at x^, and inversely 
proportional to the distance xj? and the exponent k of the attenuation of 
the source strength. 
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